UNCLASSIFIED 


AD  NUMBER 


AD480108 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  JAN  1966. 
Other  requests  shall  be  referred  to  Air 
Force  Aero  P{ropulsion  laboratory, 

Wright -Patterson  AFB,  OH  45433. 


AUTHORITY 


AFAPL  ltr,  12  Apr  1972 


THIS  PAGE  IS  UNCLASSIFIED 


480108 


AFAPL-TR-66-12 


UNCLASSIFIED 


CENTRIFUGAL  PUMP  (HIGH  PRESSURE) 
FOR  POWER  TRANSMISSIONS 


Paul  Hildebrand 
Tony  Vodopia 
Joseph  Sanders 
Lorant  Nagyszalancsy 


AIRESEARCH  TECHNICAL  REPORT  NO.  AFAPL-TR-66-12 
January  1966 


Prepared  under  Contract  No.  AF  33(657) - 1 0 1 5 1 
AiResaarch  Manufacturing  Co.Tipi.ny,  Division  of 
The  Garrett  Corporation 

for 

United  States  Air  Force 
Propulsion  and  Power  Branch 
Aerospace  Power  Division 
Wright-Patterson  Air  Force  Base,  Ohio 


UNCLASSIFIED 


NOTICES 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  what  soever;  and  the  fact  tnat  the 
Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the 
said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by 
implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 


Copies  of  this  report  should  not  be  returned  to  the  Research  and 
Technology  Division  unless  return  is  required  by  security  considerations, 
contractual  obligations,  or  notice  on  a  specific  document. 


AFAP  L-TR-66- 12 


UNCLASSIFIED 


CENTRIFUGAL  PUMP  (HIGH  PRESSURE) 
FOR  POWER  TRANSMISSIONS 


Paul  Hildebrand 
Tony  Vodopia 
Joe  Sanders 
Lorant  Nagyszalanczy 


UNCLASSIFIED 


FOREWORD 


This  report  is  the  culmination  of  a  High  Pressure,  High  Temperature 
Centrifugal  Pump  Program  conducted  by  AIResearch  Manufacturing  Company,  a 
division  of  The  Garrett  Corporation,  Los  Angeles,  California,  under  Research 
and  Technology  Division  Contract  No.  AF  33(657)- 101 51 .  The  program  was 
designated  Task  No.  8128-07  under  Project  No.  8128,  and  covered  the  period 
from  16  January  1963  through  15  December  1965. 

The  project  was  administered  for  the  Research  and  Technology  Division, 
Wright-Patterson  Air  Force  Base,  Ohio,  by  the  Air  Force  Aero  Propulsion 
Laboratory.  At  AiResearch,  this  project  was  under  the  administrative  con¬ 
trol  of  R.  L.  Schinnerer,  Chief  of  Environmental  Control.  J.  C.  Riple, 
Project  Engineer  was  responsible  for  overall  guidance  and  coordination  under 
the  direct  supervision  of  L.  R.  Woodworth,  Senior  Project  Engineer.  Major 
AiResearch  contributors  were  N.  Van  Le,  Chief  of  Aerodynamics,  L.  T.  Sladek, 
Project  Engineer,  Turbomachinery  Project,  and  J.  W.  Meermans,  Supervisor  of 
Design,  Turbomachinery  Project.  This  report  was  submitted  January  22,  1966. 


The  AiResearch  Document  Number  for  this  Program  Report  is  660095. 

Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the 
icport's  findings  or  conclusions.  It  is  published  only  for  the  exchange  and 
stimulation  of  ideas. 


Paul  W.  Montgomery,  Major,  USAF 
Chief,  Propulsion  and  Power  Branch 
Aerospace  Power  Division 
Wright-Patterson  Air  Force  Base,  Ohio 


ABSTRACT 


This  report  describes  the  analytical  study,  design,  and  experimental 
investigation  of  a  four-stage  centrifugal  pump  capable  of  30  gpm  flow  at  a 
pressure  §f  4000  psi  when  pumping  hydraulic  fluid  at  an  inlet  temperature 
of  600° F. 

The  hydrodynamic  analysis  includes  the  computer  program  format  predicting 
performance  of  elements  of  different  pump  configurations.  Test  results  of 
62  percent  efficiency,  including  all  mechanical  losses, conf i rmed  the  validity 
of  the  analysis  and  pump  design.  A  dynamic  floating  carbon  bushing  inter¬ 
stage  seal  element  was  developed  and  is  described. 
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SECTION  I 


PROBLEM  STATEMENT 


OBJECTIVES 

The  primary  objective  of  the  research  and  development  program  was 
to  advance  the  state  of  the  art  of  small,  low  capacity,  high  pressure 
(low  specific  speed)  centrifugal  pump  technology.  Specific  require¬ 
ments  were  to: 

a.  Theoretically  derive  an  impeller  and  diffuser  configuration 

applicable  to  a  pump  operating  at  45,000  rpm  delivering  30  gpm 
of  MLO  8200  hydraulic  fluid  or  equivalent  at  4000  psi  with  a 
minimum  overall  efficiency  of  65  percent.  The  fluid  inlet 
conditions  are  a  maximum  of  100  psi  and  600°F. 

Investigate  pump  characteristics,  configurations,  and  problem 
areas. 

c.  Fabricate  a  breadboard  pump  and  conduct  an  experimental 
investigation  on  the  behavior  of  the  pump  under  various 
influences. 
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SECTION  2 


SUMMARY 


ACCOMPLISHMENTS 

This  report  presents  the  results  of  a  theoretical  and  experimental 
investigation  into  the  use  of  a  centrifugal  pump  in  a  low  flow,  high 
pressure,  high  temperature  hydraulic  system. 

Hydrodynamic  Design 

The  hydrodynamic  design  studies  led  -to  the  selection  of  a  four-stage 
centrifugal  pump  with  60-degree  backward  curved  blades.  The  pump  is  designed 
to  deliver  a  flow  rate  of  30  gpm  at  a  pressure  of  4000  psi  when  pumping 
MLO  8200  hydraulic  fluid  at  a  fluid  inlet  temperature  of  600°F.  The  pump 
operates  at  a  speed  of  45,000  rpm. 

The  analytical  design  solution  was  based  on  the  calculation  results  of  a 
high  speed  computer  program  used  to  predict  the  performance  of  various  pump 
configurations  and  elements  over  a  complete  range  of  speeds.  The  calculation 
procedures  and  computer  program  format  are  included  in  this  report. 

Mechanical  Design  and  Experimental  Performance 

Centrifugal  pump  hydrodynamic  elements  were  designed  in  accordance  with 
trie  results  of  the  analytical  studies.  A  single-stage  test  rig  was  fabri¬ 
cated  to  confirm  the  hydrodynamic  design.  Results  of  the  calibration  of 
the  60-degree  backward  curved  impeller  indicated  an  impeller  efficiency  of 
approximately  80  percent  with  an  overall  stage  efficiency,  including  mechan¬ 
ical  losses,  of  65  percent. 

An  interstage  seal  suitable  for  stage  pressures  of  1000  psi  and  fluid 
temperatures  of  600°F  to  700°F  was  successfully  developed. 

A  four-stage  pump  was  designed  and  fab,  cated.  The  pump  operated  for 
a  period  of  over  twenty  hours,  of  which  six  hours  and  30  minutes  were  at  a 
temperature  of  600C'F  fluid  inlet  temperature.  The  design  provided  for  self 
lubrication  of  all  bearings.  Performance  tests  of  the  four-stage  pump 
indicated  an  overall  efficiency  of  approximately  62  percent  which  includes 
all  losses  attributable  to  shaft  seals,  interstage  seals,  and  pump  bearings. 

CONCLUSIONS  AND  RECC  MENDATI0NS 

Experimental  resuits  of  the  pump  have  confirmed  the  validity  of  the 
design  and  of  the  computer  program  pred:ctions.  The  tests  confirm  the 
capabilities  of  high-speea  centrifugal  pumps  as  hydraulic  components 
i n  a i rborne  systems. 
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The  experimental  investigation  indicates  that  the  further  development 
is  necessary  to  produce  mechanical  components  capable  of  operating  for 
extended  time  periods  in  MLO  8200  fluid  at  fluid  temperatures  of  600°F 
and  higher.  Further  testing  of  the  pump  is  recommended  in  order  to  establish 
its  mechanical  integrity  and  to  develop  mechanical  components  with  low 
parasitic  losses. 


SECTION  3 


HYDRODYNAMIC  DESIGN 


FEASIBILITY  STUDY 

At  the  initiation  of  the  high  temperature  hydraulic  purn^  development 
program,  a  feasibility  study  was  conducted.  Hydrodynamic  calculations 
were  coupled  with  centrifugal  pump  design  and  test  experience  in  an  effort 
to  assure  accurate  comprehensive  results.  The  procedure  for  t-’s  initial 
study  was  as  fol lows. 

The  centrifugal  pump  to  be  investigated  was  characterized  by  a  low 
specific  speed,  a  small  size,  and  a  requirement  fcr  high  performance. 

In  compliance  with  the  specification,  this  pump  was  to  be  designed  for 
operation  with  MLO  8200  hydraulic  fluid,  delivering  a  flow  of  30  gpm  at 
4000  psi.  Maximum  allowable  Inlet  pressure  was  100  psia  and  the  temperature 
could  be  as  high  as  600°F.  Operating  speed  was  also  to  be  taken  at  45,000 
rpm. 


This  feasibility  study  was  based  on  a  recent  AiResearch  investigation 
on  a  high  performance  small  capacity  pump.  It  was  also  conducted  with  the 
AiResearch  background  of  proven  high  performance,  sma 1 '  capacity  centri¬ 
fugal  compressors.  Latest  advancements  in  the  state  of  the  art  of  these 
compressors  has  recently  been  translated  for  application  to  centrifugal 
pumps,  which  are  indeed  old  in  theory  and  in  application. 

Efficiency  Versus  Number  of  Stages 

In  order  to  carry  out  this  parametric  study,  pump  performance  data 
were  gathered,  based  on  actual  tests  of  a  large  number  of  pumps.  This 
data  is  ^resented  In  Figure  I  .  The  attainable  efficiency  of  these  pumps 
is  shown  as  a  function  of  the  specific  speed  for  a  Reynolds  Number  larger 
than  I07.  To  account  for  the  operation  of  smaller  pumps  in  a  viscous  oil 
presently  considered,  a  correction  was  introduced  using  the  data  of  Refer¬ 
ence  I?.  By  this  method,  the  attainable  optimum  performance  of  pumps 
operating  under  various  Reynolds  Numbers  is  presented  in  Figure  2  . 

Since  the  pump  involved  in  this  study  was  small,  scaling  effects 
;  st  be  considered. 

One  approach  to  arriving  at  a  required  pump  characteristic  is  to  use 
th'  similarity  or  scaling  law  of  turbo-machines.  Derivation  of  the  per¬ 
formance  of  pumps  of  a  given  size  can  be  made  starting  from  a  proper  pro¬ 
totype. 
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Figure  2.  Estimated  Pump  Efficiency  Versus  Specific 
Speed  for  Vd^ious  Reynolds  Numbers 
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In  scaling  from  a  mode!  pump  to  a  new  requirement,  the  flow,  speed, 
and  diameter  are  given  by  similar  laws  as  follows: 


FF  “  (iiV) 

v  '  m 


Pumps  of  various  sizes,  designed  with  the  same  hydrodynamic  features 
and  same  performance  parameters  will  have  the  same  performance,  provided 
the  mechanical  design  features  and  manufacturing  tolerances  can  be  main¬ 
tained  within  the  same  proportion.  For  instance,  a  seal  clearance  of 
0.002  inches  on  a  10-inch  wheel  must  be  maintained  at  0.001  inch  to  achieve 
good  performance  with  a  5-inch  wheel. 

If  a  5-in.  impeller  has  the  same  Reynolds  number,  say  I06,  as  a  10- 
inch  impeller  manufactured  with  a  surface  finish  of  RMS  32  (root  mean 
square  of  32  microinches),  the  small  impeller  must  be  made  with  a  finish 
of  RMS  16  to  duplicate  the  performance  of  the  large  pump. 

Typical  penalty  in  performance  is  presented  in  Figure  3  .  This  curve 
has  been  established  as  a  result  of  comparing  actual  test  data  with  theo¬ 
retical  performance  for  many  different  impellers. 

The  pump  must  perform  at  a  wide  range  of  density  and  viscosity  con¬ 
ditions  because  of  the  wide  operational  temperature  spectrum.  It  will 
operate  with  a  low  Reynolds  number  (200,000)  at  low  temperature  and  still 
should  deliver  a  larger  head  with  lower  dens i ty" f lu id  at  600°F  operation. 
Because  of  these  considerations,  the  feasibility  study  was  conducted  with 
the  design  point  viscosity  at  80°F  but  the  density  at  600°F. 

On  this  basis,  the  performance  of  this  pump  was  computed  while  oper¬ 
ating  at  45,000  rpm  utilizing  different  numbers  of  stages.  Results  are 
presented  in  Figure  4  ,  indicating  that  the  pump  should  have  at  least  3 
stages  to  produce  the  desired  performance.  To  have  a  safety  margin,  a 
4-stage  design  is  preferable. 

Influence  of  the  Operating  Speed 

Although  the  problem  statement  specifies  a  speed  of  45,000  rpm,  a 
short  investigation  was  performed  to  determine  performance  at  different 
operating  speeds.  The  inlet  and  discharge  diameters  were  also  calculated. 
Figure  5  shows  the  change  of  the  inlet  diameter  at  different  operational 
speeds.  The  inlet  diameter  is  not  a  function  of  the  number  of  stages, 
but  it  does  depend  on  the  inlet  vane  angle  which  must  be  chosen  to  meet 
the  suction  specific  speed  requirements  to  assure  a  cavitation  free  operation. 
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Figure  3.  Effect  of  Impeller  Size  on  Pu^p  Efficiency 
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In  Figure  6,  the  impeller  tip  diameters  are  presented  at  different 
number  cf  stages  as  a  function  of  the  speed. 

The  calculated  efficiencies  are  shown  in  Figure  7  for  different 
numbers  of  stages.  It  can  be  seen  that  the  best  performance  is  near 
the  specified  45,000  rpm  at  a  higher  number  of  stages. 

Effect  of  Blade  Angle  on  Performance 

The  performance  of  the  low  specific  speed  machines  is  generally  lower 
than  that  of  the  higher  specific  speed  ones.  This  is  due  to  the  increased 
friction  of  the  impeller.  At  such  low  Reynolds  number  operation  like  the 
contemplated  desion,  an  additional  problem  arises,  namely  the  decreasing 
efficiency  of  the  diffuser  because  Of  the  increasing  boundary  layer  and  in¬ 
creasing  friction  factor  and  separation  effect.  To  be  able  to  choose  the 
highest  performance  design  for  this  application,  a  one-dimensional  inves¬ 
tigation  was  made  by  varying  the  discharge  vane  angle  of  the  impeller.  A 
4-stage  unit  was  chosen  for  this  investigation  with  an  operating  sp«_?d  of 
45,000  rpm  at  an  inlet  temperature  of  80CF.  The  inlet  diameter  was  kept 
constant;  the  tip  diameter  chang-d  with  the  slip  factor.  The  slip  factor 
was  calculated  according  to  B.  Eck' s  equation  (Reference  5)  which  is  a  modi¬ 
fied  form  of  Stodola1 s  formula.  It  is  in  good  correlation  with  8useman'  s 
and  Weinig's  investigations.  Furthermore,  it  allows  the  investigation  of 
the  forward  curved  vaned  impellers  in  general.  Figures  8  through  !C  show 
the  calculated  h  -ad  coefficient  change  with  the  vane  angle  at  different 
number  of  vanes  and  at  different  flow  factors  as  parameters. 

These  curves  show  in  general  that  the  head  factor  (and  the  head)  in¬ 
creases  going  to  lower  flows  below  90°  vanes  but  decreases  at  bigger  angles 
(forward  blades).  This  means  that  the  shuteff  head  is  higher  using  lower 
vane  angles,  so  tf  e  operating  characteristics  become  more  stable.  The 
curves  also  show  a  definite  increase  in  the  head  factor  using  a  greater 
number  of  vanes.  This  one-dimensional  calculation  did  not  include  an  in¬ 
ternal  loss  investigation  inside  the  impeller,  so  at  the  time,  it  was  not 
known  what  the  numerical  effect  of  the  increasing  number  of  vanes  would  be 
in  percentage  of  the  total  efficiency.  However,  the  higher  number  of  vanes 
give  little  improvement  in  efficiency  even  when  the  same  impeller  loss 
factor  was  used  in  every  case  (see  Figure  II).  The  most  promis'nn  design 
was  approximately  8  vanes  for  smaller  than  90°  vane  angles  and  16  (£  blades 
and  8  splitter  vanes)  for  bigger  angles.  Recent  small  pump  experience 
shows  that  an  8-vaned  90°  vane  angle  impeller  had  no  measurable  efficiency 
decrease  due  to  the  increased  number  of  vanes. 

Figure  12  presents  the  calculated  efficiencies  of  an  8-vaned  impeller 
versus  vane  angle  and  the  flow  factor  as  a  parameter.  Obviously,  the 
efficiency  of  th  pump  does  not  change  too  much  at  relatively  low  flow 
actors.  The  calculated  best  performance  is  with  a  0.05  flow  factor,  but 
operation  at  such  a  low  flow  requires  very  small  flow  angles  when  entering 
the  diffuser,  which  results  in  a  long  flow  path  and  fast  boundary  layer 
growth  in  the  diffuser  (Figure  13).  This  effect  was  not  included  in  the 
present  one-dimensional  calculations,  but  it  was  later  investigated  during 
the  actual  three-dimensional  computer  design  calculations.  The  difference 
in  efficiencies  between  0.05  and  0.10  flow  factors  is  only  about  5  percent. 
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Figure  6.  Discharge  Diameter  Versus 
Speed  and  Number  of  Stages 
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Figure  9  .  Head  Factor  Versus  Vane 
Angle  -  8  Vanes 
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Fiqure  10  .  Head  Factor  Versus  Vane 
Angle  -  16  Vanes 
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12.  Efficiency  Versus  Vane  Angit 
At  Different  Flow  Factors 


Figure  13 


Diffuser  Flow  Angie  Versus  Flow 


Pi ffusing  System 


The  success  of  the  program  depended  upon  achieving  a  high  perform¬ 
ance  diffusing  system.  From  this  feasibility  study,  it  was  found  that 
the  problem  was  to  design  a  diffuser  of  small  heiynt  (0.050-0.100  in.) 
capable  of  handling  a  low  Reynolds  number  of  the  order  of  4,000-16,000 
(based  on  the  height  of  the  vaneless  space).  In  addition,  this  diffuser 
had  to  readily  allow  the  transfer  of  the  flow  from  one  stage  to  the  other. 

From  the  results  of  various  diffuser  investigations,  the  following 
conceptual  designs  appeared  premising; 

1 .  Diffusing  System 

This  system  included  a  very  short  section  of  vaneless  space  (0.030- 
0.040  in.)  between  the  impeller  tips  and  the  guide  vanes  to  avoid  mixing 
losses  due  to  the  wake  of  the  blade.  A  set  of  4  to  8  discrete  diffusers 
is  shown  in  Figure  14.  These  diffusing  passages  are  designed  with  the 
proper  rate  of  diffusion  to  handle  the  thick  secondary  layer  as  found  at 
the  inlet.  In  addition,  three-dimensional  diffusion,  i . e. ,  diffusion  in 
the  axial  direction,  must  be  considered  to  obtain  a  favorable  cross  sec¬ 
tion,  in  order  to  accomodate  secondary  flow  development.  Every  diffusing 
channel  is  followed  by  a  crossover  tube  designed  to  avoid  secondary  losses. 
This  crossover  tube  leads  to  the  inlet  of  the  next  stage. 

2.  Transition  via  a  Mushroom  Diffuser  and  Plenum 


As  an  alternate  to  the  concept  of  the  crossover  ducts,  the  flow  at 
the  exit  of  the  diffusing  channel,  already  at  the  initial  stage  of  separ¬ 
ation  in  a  we  I  I -des i gned  channel,  is  allowed  to  diffuse  in  a  mushroom 
type  of  diffuser,  Figure  15,  and  then  dumped  into  the  plenum  with  very 
small  velocity.  As  a  result,  there  exist  no  swirls  in  the  plenum  that 
would  be  accelerated  at  the  inlet  to  the  next  stage,  an  effect  that  intro¬ 
duces  considerable  development  in  matching  a  multi-stage  pump. 

3.  Guide  Vanes  -  Annulus  Space  and  Returning  Vanes 

The  diffusing  system  includes  these  elements  (Figure  16)  in  a  stan¬ 
dard  approach  to  a  multi-stage  pump  design.  It  does  not  hold  too  much 
promise  for  c  low  Reynolds  number  or  small  height  application  since  the 
opportunities  for  mismatching  these  elements  are  considerable. 

4 .  Guide  Vanes  with  Boundary  Layer  Control 

To  handle  the  low  Reynolds  number  and  distorted  flow  of  the  diffuser, 
boundary  layer  control  can  be  applied.  Promising  methods  are: 

a.  Vortex  Generator  -  The  generation  of  vortices  has  been  found 
to  be  helpful  in  the  diffusion  process.  Concep tua 1  I y ,  a  small  vane  could 
be  inserted  in  the  middle  of  the  blade  channel  for  this  purpose  Figure 
17. 
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figure  14  .  Diffusing  Channels 
a'  Crossover  Ducts 


Figure  17.  Vortex  Generator 

and  Cusped  Diffuser 
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b.  The  Cusped  Diffuser  -  Preliminary  investigation  by  Refer¬ 
ence  23  shows  that  use  of  a  cusped  design  is  most  helpful  in  handling 
tht  diffusion  of  distorted  flow.  Physically,  such  a  diffuser  is  also 
shown  in  Figure  17. 

Results  of  the  Feasibility  Study 

As  a  conclusion  of  this  one-dimensional  investigation,  it  appeared 
that  a  60°  vaned  impeller  (backward  blading)  with  8  blades,  and  a  120° 
vaned  impeller  ( forward  blading)  with  8  blades  and  8  splitter  vanes  are 
the  two  most  promising  solutions.  The  backward  curved  concept  offered 
a  more  stable  head-capacity  curve.  The  forward  curved  impeller  has  a 
smaller  disc,  and  probably  better  impeller  performance.  It  has  a  lower 
degree  of  reaction  and,  therefore,  emphasizes  the  achievement  of  a  high 
performance  diffuser. 

Upon  completion  of  the  preliminary  feasibility  study,  a  generalized 
computer  program,  relating  to  all  centrifugal  pumps,  was  established. 

The  conceptual  theory  and  required  input  data  for  this  program  is  re¬ 
viewed  below  in  the  section  that  follows. 
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COMPUTER  PERFORMANCE  CALCULATION  FOR  CENTRIFUGAL  PUMPS 


This  section  presents  the  basis  for  calculating  the  performance  of 
centrifugal  pumps.  The  calculation  is  carried  out  with  a  computer  program 
that  handles  the  problem  of  performance  at  design  point  as  well  as  that  of 
of f-des i on  operation. 

DESIGN  POINT  PERFORMANCE 

Calcul  at. r  ,vi  of  De s  ign  Geomet ry  and  P e rformance  at  Design  Point 

Th>  irst  part  of  the  calculation  investigated  the  performance  of 
differed  design  combinations  capable  of  meeting  a  given  problem  statement 
as  expressed  by  the  flow  and  pressure  rise  required  by  the  fluid.  The  geo¬ 
metric  ..  ii  tens  ions ,  as  well  as  the  losses  and  efficiencies  of  these  combina¬ 
tions,  at  i  presented  to  the  designer,  who  can  then  choose  the  one  most  suit¬ 
able  for  lis  purposes. 

'■he  program  is  based  on  one-dimensional  calculation  of  the  flow  along 
the  mean  streamline,  including  friction,  diffusion,  loading,  and  boundary- 
layer-bc i Idup  factors.  The  equations  contain  experimental  factors,  which  can 
be  adjusted  whenever  the  state  of  the  art  requires  it. 

The  basic  geometric  design  form  of  the  pump  is  that  of  the  Francis  type, 
con:  sting  of  a  centrifugal  impeller  having  an  integral  mixed  flow  inducer 
portion  and  a  radially  arranged  diffuser  followed  by  a  collector  torus  or 
scroll.  This  design  arrangement,  the  most  efficient  combination,  is  shown 
in  Figure  18  . 

Within  this  general  layout,  there  are  many  possible  detail  designs.  The 
present  program  calculates  four  different  inducer  variations  and  four  differ¬ 
ent  diffuser  designs. 

I .  Input  Constants 

The  general  design  point  performance  program  requires  the  following  input 
values^  which  will  be  kept  .  in var i an t  during  the  cal cu I  at ion§ : 

Flow,  V  -  ftVsec 
o 

Pressure  rise,  &P  -  psi 
o 

Density,  YQ  -  Ib/cu  ft 
Kinematic  viscosity,  -  ft*/sec 
Inducer  design  case  (Case  I  to  4) 

Diffuser  design  case  (Case  I  to  4) 

Impeller  loss  factors  (^  t  ,  ?j,  £*  ) 
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The  necessary  Inputs  tor  the  inducer  and  diffuser  design  cases  are 
different,  depending  on  the  case. 

2 .  Input  Variables 

The  performance  of  a  pump  impeller  depends  on  (l)  its  geometric  configu¬ 
ration,  (2)  the  discharge  vane  angle,  (3)  the  slip,  which  is  primarily  a 
function  of  the  number  of  vanes,  and  (4)  the  discharge  flow  angle  or  flow 
factor.  After  these  factors  have  been  varied,  the  combination  that  performs 
best  can  be  selected.  In  many  cases,  the  optimum  running  speed  is  still  to 
be  established,  so  that  the  speed  is  also  variable. 


The  variable  inputs  are  then 


Impeller  tip  vane  angle,  B 
Impeller  slip  factor,  u  - 


Impeller  flow  factor,  \  = 


1  »  deg 


C*M 


u* 


Speed,  N  -  rpm 

3.  Basic  Relations 


The  basic  relations  used  to  establish  the  performance  and  geometrical 
characteristics  of  the  pumps  are  derived  in  the  following. 

Pump  Inducer 

The  size  of  the  pump  inducer  is  governed  essentially  by  the  availability 
of  the  suction  head,  A  high  suction  head  permits  the  flow  to  accelerate 
through  the  inducer  without  danger  of  cavitation.  The  inducer  size  is  then 
reduced  due  to  high-speed  flow.  With  the  condition  of  low  NPSH  (net  positive 
suction  head)  the  danger  of  cavitation  is  immediate,  the  flow  speed  must  be 
kept  low,  and  the  inducer  must  be  much  larger.  Thus,  the  following  four 
al ternates  are  considered  for  the  sizing  of  the  Inducer: 

1.  The  NPSH  is  given,  and  no  cavitation  at  all  is  allowed. 

2.  Inlet  vane  angle  is  given,  and  no  cavitation  is  allowed. 

3.  Inlet  vane  angle  is  given  as  60  degrees;  this  case  calculates 
low-pressure  compressors,  where  compressibility  effect  is 
neg  t i g i b I e  . 

4.  There  is  cavitation  at  the  inducer  and  inducer  dimensions  are 
precal cu I ated  . 
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Cose  I  -  NPSH  Given  and  No  Cavitation 


The  program  calculates  the  Inlet  dimensions  on  the  basis  of  the  assumption 
that  half  of  the  NPSH  can  be  used  to  produce  the  Inlet  velocity  head,  the  other 
half  of  the  head  assuring  a  cav I  tat ion-f rae  Impact  on  the  blades.  This 
assumption,  together  with  continuity  requirements,  yields  the  Inlet  diameter 
as 


576  v 

D  *  —  .  —  0  —  +  D  *  (inches) 

1  ,i  ./  32.1  75  NPSH  H 

Inlet  flow  angle  I s 


cot  P  | 


720  J  32.175  NPSH 

it  N  D. 


The  velocities  become 

C,  ■->  /  32.175  NPSH  (fps) 

I  NVD  1 

W,  =.  V  32.175  NPSH  3  (fpO 

Case  2  -  Inlet  Vane  Angle  Is  Chosen  -  No  Cavitation 

The  conventional  pump  Inlet  can  be  designed  with  this  case,  selecting 
67.5°  for  blade  angle;  other  values  can  be  used  just  as  well. 

The  Inlet  diameter  Is 

f  44,355  V 

Di  '  •  nTj'covV,  (,nches) 

where  Kj  represents  the  restriction  due  to  the  bub 


The  velocities  are  given  as: 
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The  velocity-head  necessary  to  accelerate  the  fluid  Into  the  Inlet  Is 


nr  an  ' 1  1  > 

o  64.35 

The  actual  NPSH  must  be  higher  than  this  value. 

Case  3  -  60-Degree  J/ane  Angle  Inducer 

This  case  Is  the  same  as  Case  2,  except  that  the  vane  angle  Is  fixed  at 
60  degrees,  because  this  value  gives  the  minimum  Inlet  relative  velocity  for 
Incompressible  flow.  It  must  be  noted  that  the  cavitation  performance  of 
such  an  Inlet  Is  not  too  good. 

Case  4  -  Cav I  tat Ing  Inducer 

In  a  cavl fating  inducer,  two-phase  flow  Is  encountered  at  the  Inlet. 

If  the  flow  Is  assumed  to  be  homogeneous,  the  density  of  the  mixture  can  be 
approximated  as 


Y 

-2-  i-  I 

Y  L 


»  A  Ah 


where  y 


density  of  the  saturated  liquid 


y  ~  density  of  the  mixture 

Ah  depression  head,  Inlet,  below  the  saturated  state  (TSH) 

o 

A  and  k  are  constant  charac ter  1st Ics  of  the  pumping  fluid 
For  LH|,  the  constants  are 

A  0.02014  k  0.4 

This  flow  function,  together  with  the  diameters  of  the  hub  and  tip  of  the 
inducer,  Is  used  to  calculate  the  velocity  and  angle  of  the  flow. 


144  V  (  C  *  > 

- 5  I  +  A  7-—  -  NPSH 

\ 64 . 3  5  / 


0 .  Nn 

77<r  <f|”' 


sin  P 


(fps) 
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The  flow  angle  is 


The  flow  factor  €j  represents  the  resistance  of  the  Inlet.  A  value  of 
0.95  can  be  chosen  if  the  Inlet  Is  an  undisturbed  coaxial  design. 

Losses  at  the  Inlet 

In  the  case  of  a  shockless  Inlet,  which  is  usual  In  normal  pump  design, 
the  program  does  not  separate  the  inlet  loss  from  the  general  impeller  losses. 
The  cavltating  Inducer  represents  another  case,  because  its  Incidence  is 
designed  to  overcome,  even  at  normal  operating  conditions,  the  effect  of  the 
very  low  angle  and  cavitation  in  the  inducer  throat.  This  operation  results 
in  additional  losses,  which  will  be  calculated  later  in  the  part-load  perform¬ 
ance  section. 

Losses  and  Design  Characteristics  of  the  Impeller 

The  performance  of  a  pump  as  measured  Ky  the  pressure  rise  is  often 
expressed  In  terms  of  energy  In  foot-pounds  per  pound  of  fluid  flow  or  in 
head  of  fluid  In  feet.  They  are  related  together  as 

AP 

H  =  144  — 2.  (ft) 

°  Yo 

To  deliver  this  energy,  the  pump  impeller  must  develop  a  head,  H 

through  a  change  in  momentum  imparted  to  the  flow  stream.  The  hydraulic 
efficiency  of  the  pump,  as  related  to  the  impeller  and  pump  is 

H 

o 

^HYD  M  . 

t  n 

To  start  the  calculation  in  the  program,  a  first  approximate  value  of 

n  of  0.9  Is  assumed.  Iteration  is  then  carried  out  using  the  proper  value 

HYu 

for  nnYO' 

I  .  dimeter 


are 
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For  convenience  in  performance  calculation,  dimensionless  parameters 
introduced  by  normalizing  all  values  of  head  and  losses  with  respect  to 

.  Thus,  the  Impeller  head  is  associated  with  the  head  rrmf f ic ient  defined 
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The  average  diameter  of  the  impeller  passage  is  calculated  as 


The  loss  coefficient  associated  wi tf  this  diffusion  is  then 


AqDIFF 


£3  Is  an  experimental  factor;  for  average  calculations,  0.1  can  be 
assumed. 


In  the  experimental  Investigation  of  the  interaction  between  the  Irpeller 
and  diffuser,  it  has  been  found  that  some  of  the  flow  that  has  passed  the 
impeller  does  not  enter  the  diffuser,  but  returns  into  the  impeller,  giving 
rise  to  a  considerable  energy  loss.  This  phenomenon  has  nc  been  fully 
investigated  (References  il  and  15).  At  present,  a  rec i rculal icn  loss  param¬ 
eter  Is  used  to  account  for  this  source  of  energy  dissipation. 
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The  best  value  of  ^4  factor  has  not  yet  been  determined  for  all  cases, 
but  the  calculations  of  shrouded  impeller  pumps  show  that  it  Is  close  to  zero. 
In  the  case  o*  an  unshrouded  Impeller,  we  arbitrarily  assume  the  value  0.05. 
The  foregoing  are  all  of  the  internal  losses  of  the  impeller. 


The  Interna!  static  efficiency  of  the  impeller  Is  calculated  as 
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D I f fuser 

An  Infinite  variety  of  diffuser  designs  Is  possible.  This  program 
Incorporates  four  different  solutions,  such  as 

1.  Conical  or  parallel-sided  slraight-vened  diffuser 

2.  Huitlscroll  circular  cross-section  diffuser 
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3.  Collector  scroll  (or  torus) 

4.  Vaneless  diffuser 

T he  program  calculates  the  transition  portion  from  the  impeller  Into  the 
diffuser  with  an  increased  friction  value.  This  Increase  of  loss  is  due  to 
she  mixing  of  the  uneven  discharge  from  the  impeller.  The  boundary- layer 
buildup  is  also  estimated  along  this  transition  and  the  recovery  of  the 
following  diffuser  portion  is  estimated  as  a  function  of  the  calculated  dis¬ 
pirit,  entenl  thickness  of  the  boundary  layer  (Reference  16). 

Case  I  -  Straight  Wall  Diffuser  (Figure  2Q>  ) 

The  straight-wall  diffuser  of  the  geometry  shown  in  Figure  20  is  a  design 
concept  that  lends  itself  readily  to  hardware  fabrication. 

For  calculation  purposes*  it  can  be  considered  to  be  made  up  of  a  scroll 
element  followed  by  a  straight  diffusing  section.  The  boundaries  between 
these  elements  appear  as  a  throat.  For  high  performance,  diffusers  are 
generally  designed  with  a  square  throat  section  so  that 

hTHR0AT  ,  t 
Cn  =  Z -  =  I  '  5 

u  dthroat 


The  average  velocity  ratio  in  the  scroll  is 


^  +  2  CD  tan  r 

when.  I-  „  1.04  +  0.69333  Co  B,  (^7^7") 


and  F(h)  is  a  function  characteristic  of  the  geometry  given  by 

^(2  +  tan  r>+?  loge  (j  - y 


F(h)  = 


J  **  +  61 


+  fi 


-  (_L_  ,  ftog 

Vtan  T  2-  3e 


B j (  Cptan  T  +  2)  +  0.02  tan  Y 


2B  s  +  6.02  tan  V 


Another  parameter  characterizing  the  flow  in  the  diffuser  is  the  Reynolds 
number 


Re 
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28, 


A1 T 
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With  D  for  turbulent  flow,  the  local  friction  factor  is 


i  - 


0.0462 
Re.  °-2 

bj 


This  local  friction  factor  can  be  used  to  a  Reynolds  number  lower  limit  of 
I04. 


For  Re,  lower  than  <0*,  an  adjustment  has  to  be  made,  since  the  computer 
b* 

will  extrapolate  (M)  for  calculation  purposes 


The  actual  vane  number  of  the  diffuser  as  given  later  is 

C. 


■i  ~ 


2tt 

CqB2  (2  +  CQtan  y) 


(20 


Practical  reasons  set  the  number  of  diffuser  vanes  between  6  and  16. 
Therefore,  the  program  checks  the  value  and,  if  it  is  different,  adjusts  the 
value  and  recalculates. 


This  geometry  represents  the  natural  deceleration  in  the  scroll, 
gives  the  loss 


It 


Aq 


SC 


j_ai  u,  u,  f,  +  rcuy  ill 

2<1  TT  L!  V  o’  J 


where  length  of  the  flow  path  At*  Is 


At  2  =  D  2 


"~0. 02  J\l  +  qi  B?  ~)  , 

* - - — +  7 *  •  04  +  T"  J  (lnches) 


and  the  hydraulic  diameter 


4  B.  C„  (2  *  C0  tanv) 

— - -  ^  mches) 


dH*  =  °2  2  +  2  Cn  +  Cn  tan  y 


The  loss  in  the  following  straight  diffuser  is  then  expressed  as 


X2  +  q2 

2q 


QrJ  $  r< 

Mj  L 


Cs  +  Cd  fj 


At  1 
dH* 


The  loss 
Figure  30.  £5 


factors  were 
=  0-2  and 


chosen  according  to  the  test  results  presented  in 
5s  -  1.5  give  good  approximations. 
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The  overall  diffuser  loss  Is  then 


IW1DIFF 

%  =  ruTVg  “  *q$c  +  ^D2 


Case  2  -  Hu  1 1 1  sc  ro 1 1  Diffuser  with  Circular  Diffusers  (Figure  22) 

As  the  figure  shows,  the  diffuser  design  consists  of  three  portions  -  a 
vaneless  space,  a  scroll  and  a  straight  circular  diffuser. 

As  presented  in  Reference  15,  the  loss  in  the  vaneless  space  is  equal  to 


where 


Aq  2-j  ■■= 


2  cos  ) 


(\  ,  cos;g;  "N  ft*  s 
V  1.0816  cosset*  J  2q 


0.52  f; 
B  i 


cos  a3  = 


(  I  +  cos  dj) 2  e  -  s in2a2 
( I  +  cos  ctj)  2  e^  +  sin2a* 


cot  d| 

q 

The  vaneless  space  is  assumed  to  have  a  diameter 
03  I .04  02 
The  width  ratio  is  given  as 


In  the  scroll  portion  of  the  diffuser,  the  diameter  ratio  is  given  by 


4  Dj  cot  dj  D< 


Cfc-M  ■  1 


cot  a3 


CSj- 0 


The  — 1 ratio  .s  calculated  from  this  equation. 


The  number  of  scrolls  is  an  important  design  parameter.  The  single-volute 
collector  design  corresponds  to  Z2  I .  A  higher  number  for  Z2  is  considered. 
In  the  present  program,  Z?  is  set  Z2  ~  4  for  a  typical  calculation. 
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The  loss  coefficient  in  the  scroll  is  then 


An  =  Hit  CtL  +  q2^  ^2  COS  Ot{'^  f\  +  Dj/Dj  ^ 

q5_<  2  V  2q  y  vDjcos  a3y  Vd4/d3  +  2-Bj  -  V 

\ 

For  the  straight  circular  diffusing  element,  the  minimum  possible  loss 
achieved  with  the  best  recovery  yields  a  loss  coefficient  of 


I  + 


Bj  COS2  Oj 
(04/Dj"-  0 


1 


A2  +  qA  cos  0t2\2 

cos2  of1 

[\  .  C ?  (o4/o3  +  1)  "1 

s 

^  2q  y  V03  cos  ot)^ 

Or-')1 

[  Rebj°-2(04/D3  -  l)j 

The 

boundary  layer  buildup 

factors  are 

C6  =  0.2  and  C,  = 

0  0314 

The 

total  diffuser  loss  in 

the  mu i t i score 

design  is  then 

A  =  Aq 2 -  3  +  Aq3- 4  +  Aq4-3 

qo 

Case  5  -  Single  Scroll  and  Diffuser 
This  case  is  similar  to  Case  2,  except  that  Z2  =  I. 

Case  4  -  Vaneless  Diffuser 


Aqj-j 


2 


cos 


__6 _ 

('oil  *  OfJ^ 
v  2  y 


( \  +  co s*  ai  \2  +  g2 

^  cos2  arjDj*  '  2q 


where 


3 f 2  (D3/D2  -  I) 


cos  a3 


(  I  +  cos  Otj)2  -  sin2  g; 
( I  +  cos  cr2 ) 2  e6  +  sin2  o2 


\ 

cot  a?  =  - 


After  the  long  va^eless  space,  a  quite  inefficient  additional  r<  .every 
is  possible  in  the  collector. 
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Here  a  fraction  of  the  kinetic  energy  of  the  velocity  leaving  the  inpe) ler 
Is  assumed  to  be  lost. 

Aqs- 1 ~  0.66  J  -  C9S 

s  *  \  Cos  aj  2q 

The  total  diffuser  loss  in  the  vaneless  diffuser  is  then 

%  =  Ad..-4  *  Aqt-, 

Hydraulic  Performance  of  the  Pump 

The  hydraulic  performance  measures  directly  the  work  which  produces  the 
pump  pressure  This  means  that  the  head  output  of  the  pump  must  be  equal  to 
the  required  value  after  all  the  internal  'esses  are  deducted 


The  hydraulic  efficiency  of  t he  pump  is  then 
I 

'TH 


HYD 


H 

_£ 

H, 


Sl  ~  Aqp 
1  4  ^SF  4  NC 


The  hydraulic  efficiency  is  now  calculated  from  the  above  equation  yith 
various  loss  parameters  which  are  obtained  by  starting  the  computation  with 
an  assumed  value  for  T)  ^  The  value  from  the  above  equation  must  be  within 

2  percent  of  the  assumed  value,  or  otherwise  the  computer  will  proceed  to 
I terate. 

Overall  Efficiency  and  Power  Input  of  the  Pump 

Besides  the  loss  in  energy  in  the  impeller  area  diffuser  passages, 
externa!  losses  must  be  considered.  These  are  due  t-  friction  from  the  hub 
disc  which  supports  the  impeller  blades  and  the  leakage  through  the  seal  clearance. 

The  disk  friction  loss  is  often  given  as 


6q 


FR 


Lu 


I  Q6  v 


1 0  Bj  \  Re j 


) 


The  experimental  factor  0.285  is  to  be  assumed  for  the  case  oc  a 

shrouded  wheel  and  C  366  for  an  unshrouded  wheel. 


A  q, 


Ll-£Q& . V 

e. 


f 


2g  -  X1 


2q 


The  factors  are 
an  unshrouded  one. 


0  for  a  shrouded  impeller  and  ??  0.5  or  0  for 


The  overall  efficiency  of  the  pui.p,  from  which  the  hydrodynamic  power 
input  is  obtained,  is  then 


H 

Out 

H. 


.1*1*1- 


Aq 


FR 
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PART-LOAD  PERFORMANCE  OF  THE  PUMP 


The  program  calculates  the  estimated  friction,  diffusion  and  incidence 
losres  of  the  design  described  by  the  design-point  performance  program  for 
variable  flow  and  speed  values.  Provision  is  also  made  for  changing  the 
fluid  density  and  viscosity;  the  geometric  design  is  fixed. 


Input  Variables 

In  addition  to  the  design  point  input,  additional  data  are  necessar 
for  the  part -load  program 

v  -  i he  viscosity  of  the  fluid,  which  can  be  different  from 

the  des i gn  v i scos i ty 

v  -  the  density  of  the  fluid,  this  can  be  also  different 


=  flow  ratio  in  steps,  defined  as  M 
constant  speed 


o  design  flow 


speed  i at i o  in  s teps 


design  speed 

NPSH  net  positive  suction  head;  can  be  different  from  the  design 

va  I  ue 

A  e  d  K  cavitation  expansion  constants 

ai ,aj,aj  incidence  loss  constants  of  the  inducer 

•it  cavitation  loss  constant 

c 

y3  d'ffuser  incidence  loss  constant 

!  1  I e r  Performance 


The  impel le'  inlet  flow  is  equal  to  the  nominal  flow  plus  the  increase 
due  to  leakage.  Expressed  by  the  flow  factor,  it  is 

H  M  *  Aq 
o  l 


The  leakage  loss  Is  taken  from  the  original  design  po;nt  calculation 
as  a  first  approx imat ion 
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With  the  change  in  <  I uw  rate  and  speed,  the  flow  factor  at  the  impeller 
changes  so  that  the  head  coefficient  is  now  calculated 

q  u  (  I  -  \  M  t  an  0 ) 

o  o 


u  = 


r  n  cos  TT 


1  *  -  JT-frrcTTf,  0  08  N 


v  c.  J 


wi  th 


/  h  2  \  2  +  ( 

Jo  o 

V  . v7*  — 


-  q)1 


T~TT 

o  o 


where  0 


- —  =  v2  cot  B. 

Uj  r,  10 


At  off-design,  an  incidence  is  encountered  at  the  inducer  inlet  so  that 
the  incidence  loss  is  given  as 


&q. 


inc 


M  , 
W  57 


where  is  an  exper imenta !  constant.  According  to  the  test  values  on  axial 
compressors,  *2  is  a  parabolic  function  of  the  incidence  angle: 

♦j  =  a,  *  aj  fi  8.  *  aj  4  B.! 


A8.  - 

f 


B.  -  P. 

I  i  o 


(absolute  value) 


cot  B,  M  cot  0 


The  average  relative  velocity  in  the  impeller  passages  is 

u?  r 


v  ■  1  *  J'2  (H2  cot2  8. 

.  T 


2  ' 


The  Reynolds  number  calculated  with  this  velocity  inside  the  >r'pe  Mer 
passages  is 


Re 


U  ? 


dHYD.  WAV 

-itV- 

>>  '  f,  »  N 


720  vN 


d|_jVn  •  s  given  by  the  design  point  calculation.  Tne  new  f  rict  ion  ’act 

0.0462 
f|  =  rTTT 


or  becomes 
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With  this  factor,  the  friction  Is  to  be  calculated 
f|  bl\  W. 


"oF 


Si 


U 


A  V 

T 


UHYD |  Ul  q 
tl i  is  to  be  taken  from  the  design  point  calculation. 

The  diffusion  factor  through  the  impeller  Is 

q  (cos  0.  t  cos  0,) 

A  I  -  Ci 


1 1  ’  W  i-  - 

2  u7  ^  (l  "  v*}  f  2vi 


where 


W. 

i 

U2 


»  v? *  cot*  8.  M* 

IO 


The  diffusion  loss  is  then 

4V  ■  (W 

And  the  recirculation  loss  becomes 

AqRC  L  M^X- 

o  o 

The  constants  Si»  Sa,  and  S*  are  the  same  as  in  the  design  point  program. 
2.  Cavltating  Performance 

The  calculation  of  the  cavltating  Impeller  performance  is  similar  to  that 
for  the  previously  described  one,  but  since  the  density  is  not  constant  in  the 
inlet,  the  inlet  velocities  are  a  function  of  the  flew. 


The  head  factor  is 

(I  -  X  M  tan  ) 


q 

M 

Cl 

0 

o 


I A 


I 


1  2Z ,  (  I  - 


V2 


*  ♦  (I  -  q ) ' 


— ? — r~c 

V2  1  S I  M  0 


c 

U2 
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The  inducer  entering  velocity  Cj  results  from  the  design  point  calculation. 
Following  the  correlation  of  the  mixture  density  given  previously,  the  density 
change  at  the  inlet  to  the  cavitating  inducer  changes  with  speed  and  flow 


rate  as 


let  M  ,  S") 
r  \  0  V 

y(at  design) 


I  +  A  (  M  2 

V  o  2  q 


i  K 


NPSH 


I  ♦  A  -  NPSH  ) 


where  the  design  point  corresponds  to  the  value  of  I  for  the  parameters  M 

*  Cto*\  is  smaller  than  the  c 


("o'  r-) 


and  N  .  When  the  inlet  velocity  head 

N 

NPSH,  no  density  changes  take  place.  Then  5i  =  I . 

The  inducer  incidence  loss  is  similar  to  the  previous  case: 


Aq 

inc 

*  ^2 

w. 

2q 

where  i)r  2 

« 

a  1  + 

a  2  Ap  |  +  a  3  A  P  | ' 

Ap{ 

z 

h 

-  PjQ|  (absolute  value) 

S, 

M0 

cot 

P1 

n  '“r  ■ 

o 

V  2 

vjl+  S, 2  H2  0  ? 

o 


The  loss  constants  can  be  assumed  to  be 


a  |  =  0 

3  j  ~  0 

a3  =0.025 
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There  is,  however,  an  additional  loss  with  cavltating  operat lon--the 
cavitation  loss: 


^qCAV  ' 


I  + 


*  ("o'  ^  •  NPSH)]  K  ■  ' 


(v*1  4  S.W*  *) 


The  value  of  i|i  can  be  assumed  to  be  between  0.5  and  1.0. 
The  impeller  skin  friction  is 

ft  Lit  wAV2 

6qsF  ’  5’  vJTT  ^ 


0.0462 

where  f,  -= 


Re 


U 


A  V 


dHYDi  WAV 

1  "  12  vJ> 


(i  +  Cij)(s,2m20o!  4  v*2) 


The  diffusion  factor  is 

& 


q  (cos  0.  t  cos  0? ) 

1  ‘  Cr  1  ST - 1 - 


2U7 


(I  -  V?)  +  2v* 


W. 


ut  •  v 


U; 

With  these  values,  the  diffusion  toss  is 

/W  2  \ 


AqDL  ? J  A 


,U2'q, 


and  the  recirculation  loss  will  be 


Aq 


RC 


?4  Sh 

o  o 


p,,  |»3,  and  f»4  factors  are  the  same  as  in  the  design  point  program. 
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Diffuser  Performance  ai  Part  Loud 

At  off-dcsign  operation,  the  diffuser  vanes  or  scrolls  suffer  incidence 
loss  on  the i r  tip: 

(q  -  q  )*  '  (X  -MX)1 

o  o  o  o 


ijl j  factor  Is  also  a  function  of  the  Incidence  angle.  As  the  angle  of 
attack  Is  normally  very  small  with  the  low  flow  angles  encountered  In  pump 
design  practice,  this  value  can  be  assumed  constant  and  equal  to  unity. 

The  diffusion  loss  can  be  calculated  by  assuming  the  same  relative  loss 
percentage  of  the  available  energy  as  was  assumed  at  the  design  point. 


M  q 

Aqn  Aqn  ■ 

uhD  mDo  q 


<»„V  ‘ 


where  *s  the  diffuser  loss  at  design  point. 

Pressure  Ratio 

The  head  and  pressure  can  be  calculated  from  the  hydraulic  losses. 
The  total  head  output  is 

HTH  32 . *75 
The  hydraulic  efficiency  is 

1  "  AqINC  “  AqCAV  ’  AqDL  "  ^qRC  “  AqDIN  "  AqD 


I  '  Aq< 


The  head  of  the  pump  becomes 
Ho  T,H  HTH 


and  the  pressure  rise  is 


Y  M 

TXT  (mi) 
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Outside  Losses  and  Overall  Efficienc 


The  wheel  friction  changes  with  the  flow  and  the  Reynolds  number: 


1000  \  M  qB2 

o  o 


(\ 


12  v* 


The  leakage  changes  al«'' 
f  0.006  V2 ' 


2q  -  \  2  M  *  i  qJ 


5*  and  5*  factors  are  similar  to  those  of  the  design  point. 
The  overall  efficiency  is  calculated  with  these  losses. 

^  1  *  Hr  *  H 

The  power  consumption  of  the  wheel  is  then 


H  H  V  y 

o  o  o 


r,  550 

The  program  prints  out  the  losses,  the  efficiencies,  the  head,  the  pressure 
and  the  horsepower  values.  In  addition,  it  prints  out  a  diagram  of  the  pressure 
rise  and  the  efficiency  versus  relative  flow  withrelative  speed  as  a  parameter. 


*1 
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GEOMETRY  AND  PARAMETER  CALCULATIONS,  COMPARISON  WITH  TEST  RESULTS,  AND 
NOMENCLATURE 

Calculation  of  the  Inlet  Scroll  of  a  Vaned  Diffuser  with  Conical  Sideplate 

Figure  23  shows  the  general  arrangement  of  the  inlet  scroll  of  a  vaned 
diffuser  with  conical  sideplate.  The  numerical  calculation  of  thes  inlet  scroll 
profile  is  based  on  the  natural  deceleration  principle.  This  means  that  the 
momentum  of  the  fluid  is  decreased  only  by  the  friction  on  the  walls.  The 
geometry  of  the  diffuser  scroll  is  shown  in  Figure  23  also. 
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Since  dM  =  frc*  (d«  -  bQ)  j  coscxdx  and  dm 
Equations  (  * 4) and (15)  results  In 


P  ~~£~  dx,  combining  these  with 


d<Cur) 

(cTT 

u 


f  It 


C^b  r 
Mo  o  o 


(I 


nd 


)  d  (d) 


The  Integral  Is  then 
I 

C  r 
u 


b 

-A  ( d  -  —  log  d)  4  constant 


In  this  expression  the  logarithmic  term  represents  the  so  called  "tongue" 
effect— the  Increased  friction  due  to  the  theoretically  high  losses  con¬ 
nected  with  the  very  sma I '  diameter  in  the  beginning  of  the  scroll.  This 
Is  only  theoretically  so.  In  pract Ice,  the  channel  always  has  a  minimum 
width— the  width  of  the  vaneless  diffuser  portion  where  the  flow  can 
adjust  Itself  and  reduce  the  relative  friction.  Therefore,  this  logarithmic 
term  of  the  equation  can  be  replaced  by  a  constant  and  an  Increased 
friction  factor  given  by  experimental  comparison  with  existing  diffuser 
results. 


C  r 
uo  o 


f2"  (r  -  r  ) 
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C  r 

b  tan 

a  o 
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o 

C  r 

uo 

r  .  , 

nf 

r~ 

r  *  B 

tan  a 

(rr  -  i) 

o  ’O 
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uo 


(i  .  M  (<  .  — L.) 


The  angular  position  of  the  cross  section  given  In  degrees  of  angle  Is 


-- 


90 


B  tan  a 
o  o 


( 


Cu  V  d 

r  Xir, 

UO  '  O 


This  equation  also  can  be  expressed  as  the  number  of  diffuser  vanes 


B 


It 


o  t  an 


4/d_  '  * 
ID 
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Estimating  Straight  Diffuser  Recovery  from  Boundary  Layer  Thickness  in  the 
Throat 

Reference  7  describes  some  tests  conducted  with  straight  and  curved 
channels,  measuririj  the  maximum  pressure  recovery  of  different  diffuser 
configurations  having  variable  boundary- layer  thicknesses  in  the  diffuser 
inlet  throat.  In  the  report,  the  actual  displacement  thickness  of  the  boundary 
layer  is  given. 

The  recovery  factor  of  a  diffuser  is  defined  as  the  rat'O  of  the  pressure 
rise  across  the  diffuser  and  the  total  inlet  velocity  energy  entering  the 
di f  fuser  : 


R 


C 


The  relative  displacement  thickness  the  boundary  layer  at  the  inlet 
throat  of  the  diffuser  is  given  as  the  ratio  of  'he  integra1  average  dis¬ 
placement  thickness  and  the  width  of  the  throat,  2£*  . 

b 


The  measured  test  results  were  plotted  In  Figure  25  for  straight  channels 
and  for  curved-elliptic  channels.  The  result  shows  the  well  known  fact  that 
the  straight  channel  has  better  recovery  than  a  curved  one. 


The  measured  performance  points  of  the  straight  diffusion  channels  can 
be  connected  by  a  linear  function.  The  general  form  of  the  loss  vs  relative 
boundary  layer  thickness  will  take  the  following  form: 


C 

l 


a  i 


a? 


is* 

b2 


For  general  calculations,  it  seems  to  be  better  to  replace  the  diffuser 
width  by  the  hydraulic  diameter  of  the  diffuser  throat 


C 


l 


•  i 


♦  a, 


The  a |  and  aj  constants  are  functions  of  the  form  of  the  throat  The 
circular  cross  section  gives  low^r  losses  and  higher  recovery  than  a  Hat 
rectangular  one- 


With  the  help  of  this  approx ima te  equal  ion  of  the  pressure  loss  factor, 
it  If  possible  to  estimate  the  loss  of  the  straight  portion  of  the  diffuser. 
The  displacement  thickness  can  be  calculated  from  the  average  length  of  a 
streamline  and  the  Reynolds  number,  while  the  hydraulic  diameter  is  given  by 
the  geometry  of  the  diffuser. 
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In  the  case  of  conical  sidewalled  vaned  diffusers,  the  length  of  the 
streamline  is  bt2,  the  hydraulic  diameter  is  d^j.  The  displacement  thickness 

with  turbulent  flow  is  a  function  of  Reynolds  number  and  length  of  the  stream¬ 
line: 


6”  ip,  LI, 

The  friction  factor  in  the  diffuser  function  Is  similarly  a  f l  ;c t ion  of  the 
Reynolds  number,  so  the  displacement  thickness  also  can  be  expressed  in  the 
fol lowing  form : 


Combining  the  two  previous  equations  gives 

V  a<  +  d^ 

The  total  loss  of  diffusion  after  the  throat  is  the  product  of  the  loss 
factor  given  above  and  the  available  velocity  energy  in  the  diffuser  throat. 
The  velocity  rat'^  between  the  impeller  tip  and  the  diffuser  throat  can  be 
calculated  as  shown  in  preceding  paragraphs.  The  velocity  energy  ratio 
becomes 

CT*  q*  ♦  A*  X* 

2g  Hth  2g  l  CHj>  V 

and  the  relative  loss  in  the  straight  diffuser  is 

Hiff  2  ■  -Hr4- ( c;")  (''***-*»f*^ 


or,  with  combined  toss  factors 


I  IFF  2 


^7 7^  (?*  ‘  fl  df^) 


iriso*.  with  Test  Results 


As  w,  emphasized  earlier,  this  program  calculates  the  per f oreta.sce  of  an 
assigned  deal  p-jmp  conf  igurat  «otn  and  does  not  compute  efficiency  of  ex'sting 
pump  gcomstries-  This  fact  must  be  kept  in  mind  when  calculations  >  rf  compared 
with  test  *ei’il  IS  for  existing  machinery 


In  der  to  compaic  the  calculation  results  with  statistical  average  and 
predicted  optimum  performances,  a  number  of  data  were  gathered  on  high-Reynolds- 
number  pump  performances  and  plotted  in  Figure  I  For  comparative  purposes, 
three  points  were  calculated  by  the  present  computer  program  using  water  as 
fluid,  assuming  no  cavitation  and  a  Reynolds  number  close  to  I07.  The  specific 
speeds  of  the  calculated  samples  were  500  to  1000  and  2000.  The  results  of 
the  calculations  were  also  plotted  in  the  diagram.  The  points  fit  the  top 
efficiency  curve  quite  well,  the  program  therefore  might  be  used  to  predict 
the  performance  of  advanced  pump  designs  with  confidence 

Nomenc lature 

Y  density  (ib/cu  ft) 

Yq  average  liquid  density 

P  pressure  ps  ia  ) 

APq  pressure  rise  (psi) 

H  head  (ft/lb/lb) 

H  actual  head  (ft) 
o 

total  head  developed  by  impeller  (ft) 

r\  efficiency 

T)o  actual  efficiency 

r^YO  hydraulic  efficiency 

0  relative  angle  (degrees) 

0.  inlet  flow  angle 

0i  inlet  vane  angle 

02  discharge  flow  angle 

020  discharge  vane  angle 

NPSH  net  positive  suction  heed  ft) 

0  diameter  (inches) 

D.  inlet  eye  diameter 

D  inlet  hub  diameter  at  0. 

H  I 

D2  impeller  tip  diameter 

V  flow  (cf s ) 

V  !  i quid  f low 
o 

N  speed  (rpm) 
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C  absolute  velocity  (foe) 

C.  inlet  velocity 

C*  discharge  velocity  at  impeller  tip 
CiM  inlet  meridional  velocity 

C/u  discharge  meridional  velocity 
n 

discharge  tangential  velocity 
W  relative  velocity  (fps) 

W.  inlet  relative  velocity  at  impeiler  eye 

W?0  discharge  relative  velocity  at  impeller  tip  ~  tneoretical 
W 2  discharge  lela.i/e  velocity  at  impeller  tip  -  actual 

U  rotational  veloc1;*  ,  os) 

U.  inlet  -national  v  iocity  at  impeller  eye 

U;  discharge  rotational  velocity  at  impeller  tip 
K  blockage  factor 

K.  inlet  h jb  D'ockage  factor 

r  'lei  contraction  factor 

K  iniet  density  function  exponent 

A  inlet  dersity  function  constant  (I/ft) 

3  total  head  factor  (Cj  /u*) 

u 

p  sUp  factor  (  Cj  /Cj  ) 

u  uoo 

A.  flow  factor  (c uu/ u 2 ) 

vj  impeller  diameter  ratio  (D./Dj) 

£i  impeller  diffusion  ratio  (Wj/W.) 

Cp  d'ffuser  throat  width  ratio 

Z  number  of  varies 

Z |  of  impeller 

Zj  of  d i f fuser 
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Lt  i 
Liz 
B  2 
b2 
d 


Re 

f 

tq 


A 

V 

a 


loss  constants 
heat 

impe  I  ier  f r ict  ion 
Is  diffusion 
5»  rec i rcu !at t  on 


5g  wheel  friction 
( «  leakage 


Mf 

diffuser  boundary  layer 

?6  - 

t  y  J 

mean  streamline  Mngth  finches) 
diffuse'  scroll  streamline  length  (inches) 
impeller  width  factor  (dj/D^  } 
impeller  tip  width  (inches) 
flow  channel  diameter  (inches) 

^HYD  ^vt^rdu’'c  diameter  (inches) 

d  diffuser  scroll  hydraulic  diameter  (inches) 
n  * 

Reynolds  number 
friction  factor 
loss  (dimens ionless) 

impeller  skin  friction  loss 

Aq^  impeller  diffusion  loss 

6qRr  impeller  recirculation  loss 

disuse-  scroll  friction  loss 


Ae*oai 


straight  diffuse:  loss 
kinematic  viscosity  (ft'/sec) 


Aq(  feakaga  loss 

Aq^  disk  friction  loss 

Aq-  total  diffuser  loss 

Zq,  inducer  inlet  incidence 
lflC  loss 


diffusion  factor- impe 1 ler 
diffuser  wall  angle  (deg) 
absolute  Mew  angle  (deg) 

(1,  inlet  flow  angle  (pre  rotdt  i  on  ) 

a?  impeller  discharge  f'ow  angle 
aj  Mow  angle  and  of  vane  less  space 
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PERTINENT  HYDRODYNAMIC  DESIGN  STUDIES 


Selection  of  Design  Points 

A  four-stage  pump  conf !gurat ion  has  been  chosen  to  deliyer  the  4000  psi 
required  by  the  problem  statement.  Fluid  incompressibility  permits  identical 
contouring  of  each  stage.  Thus,  comparative  :ovest igct ion  and  test  can  be 
made  for  one  stage  only. 

Two  different  designs  reviewed  in  the  feasibility  study  were  chosen  for 
fabrication  and  test  for  this  program.  The  basis  of  selection  **.’as: 

a.  Performance  at  the  design  point 

b.  Performance  at  off-design  conditions 

c.  Size  and  weight  considerations 

Two  generalized  calculations  were  established  and  put  into  the  IBM 
computer  program.  The  first  program  calculated  the  impeller  and  diffuser 
dimensions  and  velocities,  pressures  and  losses  associated  therewith  at  the 
design  point.  The  second  program  calculated  off-design  condition  performance 
based  on  the  above  dimensions.  The  partial  load  performance  calculation  was 
repeated  using  different  density  and  viscosity  values  representing  a  second 
fluid  temperature. 


The  specified  design  conditions 
Fluid 

Design  Flow 
Design  pressure 
NPSH  (first  stage) 
Temperature  range 
Density  range 
Kinematic  viscosity  range 
Design  speed 


M10  8200 
30  gpm 

975  ps i  per' stage 
100  ps  1 
80  to  600°F 

56.7  to  45.3  lb  per  cu  ft 
3.49  x  10"*  to  0.198  x  10" 
45, 000  rpm 


ft2  per  sec 


Performance  at  Design  Point 

The  fluid  conditions  at  600°F  were  conservatively  selected  as  a  basis  for 
design  point  calculations.  The  lower  density  at  this  temperature  requires  a 
h;gher  head  pump  output  to  attain  the  specified  pressure. 


The  inlet  geometry  of  the  pump  was  maintained  constant  for  each  of  the 
designs  invest igated,  because  cavitation-free  operation  can  be  achieved  at  the 
minimum  inlet  velocity  point.  The  inlet  dimensions  are  determined  by  the 
design  flow,  the  NPSH  of  the  first  stage,  and  the  minimum  hub  diameter,  as 
dictated  by  structural  considerations  for  the  four-stage  pump.  The  hub  diameter 
used  was  0.660  inch. 

In  addition  to  the  fixed  inlet  geometry,  the  variables  fed  into  the  IBM 
program  were: 

Vane  angle  at  impeller  tip  60,  30,  0  and  -30° 

Slip  Factor  (C  _/C  )  0.6  to  0.9 

K  u2  uoo 

Impeller  tip  flow  factor  0*05  to  0.25 

Program  output  was  as  follows: 

Impeller  and  diffuser  dimensions  and  velocities 
Impeller  and  diffuser  vane  numbers 
Reynolds  numbers  and  friction  factors 

Boundary  layer  buildup  in  the  channels  (based  on  one-dimensional 
veloc it ies) 

Friction,  diffusion  and  leakage  losses 
Hydraulic  efficiency  and  pressure  rise 
Overall  efficiency 

The  results  were  plotted  and  evaluated  based  on  dimensions,  performance 
and  design  simplicity. 

Figures  25  and  26  ,  included  herein,  show  the  resultant  performance 
curves  of  the  60°  backward  curved  vaned  pump  impeller  and  the  -30°  forward 
curved  configuration,  respectively.  The  best  performance  was  obtained  using 
12  and  30  impeller  blades  for  the  respective  vane  angle  configurations.  There¬ 
fore,  these  arrangements  are  shown  in  the  figures. 

The  efficiency  shown  is  computer  calculated,  and  is  necessarily  based  on 
an  arbitrary  diffuser  design  which  is  parallel  walled  and  necessarily  has 
higher  friction  and  boundary  layer  buildup  than  a  diffuser  with  conical  side 
walls  or  increasing  cross-sectional  area.  The  reason  for  using  the  parallel- 
walled  diffuser  in  the  comparative  computer  program  was  that  it  permitted 
solution  of  the  one-dimensional  flow  equations  and  loss  functions  in  a  closed 
loop  form,  as  necessary  for  the  general  computer  program. 
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Figure  25.  Pump  Performance  at  Design  Conditions 
( IBM  Calculated) 
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Off-Design  Performance 

To  fully  ascertain  pump  performance,  an  off-design  performance  analysis 
was  also  conducted,  using  600°F  and  80°F  fluid  temperatures. 

The  design  condition  pump  geometry  was  reused  for  the  off-design 
performance  program,  which  determined  pressures  and  losses  at  various  fluid 
flow  rates  and  viscosities. 

Results  are  presented  in  Figures  27  and  28  ,  which  show  pump  head  and 
efficiency  as  functions  of  flow  for  the  60°  backward  curved  and  -30°  forward 
curved  impeller  versions.  The  room-temperature  performance  run  resulted  in 
higher  pressures  but  lowtr  efficiencies  than  for  hot  operation  due  to  increased 
density  (and  a  lower  Reynolds  number)  at  the  increased  viscosity.  This 
pressure  increase  can  be  reduced  by  decreasing  pump  speed  to  41,000  rpm.  This 
is  easily  accomplished  if  the  pump  were  turbine  driven.  In  the  present  case, 
the  dynamometer  driving  the  test  unit  can  readily  handle  this  speed  change. 

On  the  other  hand,  the  increased  viscosity  can  produce  higher  system  resist¬ 
ance,  in  which  case  the  higher  pressure  may  be  desirable. 

Relative  Merits  of  Backward  and  Forward  Bladed  Designs 

Comparison  of  the  backward  and  forward  curved  impeller  curves  revealed 
that  the  former  has  a  definite  superiority  in  design  point  efficiency  and 
better  pressure  curve  stability.  The  3.5-percent  calculated  efficiency  differ¬ 
ence  is  due  to  the  higher  calculated  losses  in  the  diffuser.  It  is  possible 
to  minimize  this  loss  with  optimum  diffuser  design,  and  boost  efficiency  to 
the  64  to  65  percent  level. 

The  main  objection  to  the  forward  curved  design  is  its  unstable  head- 
capacity  curve.  The  curve  presented  in  Figure  28  is  the  result  of  the  computer 
calculations  using  average  design  parameters.  It  is  certain  that,  with  care¬ 
ful  design  of  the  components,  the  efficiency  can  be  increased  and  the  head- 
capacity  curve  can  be  stabilized.  It  should  be  added  that  the  smaller  and  low 
friction  producing  forward  curved  pump  design  has  greater  performance  improve¬ 
ment  potential. 

Detailed  Design  Consideration' 

I .  Inlet  Channel 

To  reduce  diffusion  in  the  impeller,  the  inlet  must  be  contoured  .0  that 
the  relative  inlet  velocity  at  the  impeller  tip  is  at  a  minimum.  A  low  specific 
speed  pump  impeller  converts  a  very  high  percentage  of  the  total  energy  into 
pressure,  so  that  efficiency  improvement  is  extremely  important.  Figure  29 
shows  the  average  relative  velocity  change  for  a  shockless  inlet  at  various 
inlet  flow  angles.  The  relative  velocity  becomes  minimal  at  an  inlet  vane 
angle  of  approx imat e I y  68°. 

The  inlet  channel  is  contoured  in  the  form  of  an  elliptic  toroid  as  shown 
in  Figure  30.  This  fo^m  is  nc i i  proven  in  various  nozzle  applications  as  one 
which  produces  smooth  velocity  distribution  both  across  the  channel  and  along 
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RADIAL  DIMENSION, 


the  shroud  wall.  The  I ovj  pressure  gradient  associated  with  the  high  velocity 
at  the  impeller  inlet  face  minimizes  the  possibility  of  flow  separation  and 
cavitation.  The  slightly  elongated  form  of  the  channel  w?s  necessary  to 
accommodate  the  inlet  seal. 

2 .  Impel ler  De*  ign 

a.  60°  Backward  Curved  Impeller  (Conventional  Design)  -  The  comparative 
IBM  calculation,  presented  in  Figure  25  ,  shows  an  optimum  design  at  the 
fcl lowing  condi t ions: 


Flow  factor  (p-1—  ) 
U2 

Slip  factor  (—■  } 
U2 


=  0.  I  I 


*  0.739 


Number  of  blades  (Z) 


=  12 


Impeller  diameter  (D^)  =  2.308  inches 

Impeller  tip  width  { b£ )  =  0.040  inch 

These  parameters  provide  the  basis  for  the  actual  design.  The  inlet  and 
the  hub  diameters  were  established  as  reviewed  in  the  previous  sections. 

The  design  was  established  in  five  steps. 

a.  The  shroud  streamline  was  established  between  the  eye  diameter  and 
the  t;p  usmg  a  natural  hyperbolic  function  to  assure  a  smooth 
change  in  curvature  and  a  big  radius  of  curvature  at  the  eye.  This 
curve  can  be  expressed  by  the  equation: 


where  r  =  the  rad'u'. 


Z  =  the  axial  coordinate 


a,  b,  c,  n  =  constants 


b.  The  vane  angle  change  along  the  shroud  was  established.  This  was 
done  using  the  conformal  t ransf ormot ion  method.  The  vane  is  estab¬ 
lished  as  a  parabolic  arc  in  the  transformed  coordinate  system, 
expressed  as  a  function: 


t  an  P  = 


=  a  §  +  b 
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and 


dm  = 

rd§ 

dn  = 

dcp 

where  0 

?  and  T] 
m 

? 

r 

a  and  b 


vane  angle 

transformed  coordinates 

arc  length  along  the  shroud  line 

angular  position 

radius 

constants 


c.  The  desirable  relative  velocity  distribution  along  the  shroud  was 
assumed  and  the  average  flow  areas  and  impeller  widths  were  computed. 
This  establishes  the  hub  line. 

d.  An  electric  analog  field  plot  was  made  to  obtain  two-dimens ional, 
rotationless  potential  lines  along  the  impeller  passages. 

e.  The  three-dimensional  velocity  distribution  in  the  rotating  passages 
was  then  calculated  by  the  IBM  computer. 

Figure  31  shows  the  resultant  velocity  distribution  along  the  shroud 
line  for  the  60°  backward  curved  impeller.  The  number  of  vanes  was  increased 
from  12  to  16  because  the  three-dimensional  calculation  showed  stagnation 
along  the  pressure  side  of  the  vanes. 

The  16  vanes  result  in  very  good  pressure  and  velocity  distribution  along 
the  passage,  but  cause  more  than  permissible  blockage  at  the  inlet.  To  avoid 
this,  the  number  of  vanes  was  reduced  to  8  for  the  inlet  portion  and  16  (using 
splitter  vanes)  near  the  outlet,  or  large  diameter,  portion  as  shown  in 
Figure  32. 

The  vane  sections  were  radial  elements  from  shroud  to  hub,  thus  assuring 
easy  mach i nab i i i ty .  The  thickness  of  the  vanes  is  0.020  in.  at  the  inlet  and 
a  constant  0.030  in.  thereafter.  Slight  tapering  of  the  vanes  toward  the  hub 
is  also  utilized  due  to  structural  considerations. 
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Figure  32. 


b.  Minus  30°  Forward  Curved  Impeller  -  Figure  26  presents  the  optimum 
efficiency  at  the  following  conditions: 


M2 

U2 

=  0.116 

CU2 

U2 

=  0.986 

Number  of  blades  (Z) 

=  30 

Impeller  diameter  (D^) 

=  1.986 

i nches 

Impeller  tip  width  ( b^) 

=  0.039 

inch 

The  inlet  tip  and  hub  diameters  are  established  as  described  previously. 

The  five  design  steps  reviewed  above  were  also  utilized  in  this  design. 

A  deviation  from  the  established  design  point  was  necessary.  The  number 
of  vanes  was  reduced  from  30  to  24  because  of  inlet  blockage  considerations- 
The  maximum  number  of  vanes  of  0.02C  inch  thickness  is  12  at  the  inlet  edge 
without  creating  a  blockage  problem.  The  resultant  reduction  in  output  pres¬ 
sure  is  calculated  to  be  less  than  0.5  percent.  Incorporation  of  intermediate 
(or  splitter)  vanes  in  the  large  diameter  portion  results  in  the  compromise 
value  of  24  vanes. 

Figure  32  presents  the  velocity  distribution  along  the  shroud  for  this 
case.  The  loading  is  higher  than  for  the  backward  curved  impeller  case,  but 
it  is  still  acceptable  and  stagnation  has  Keen  avoided. 

The  vane  sections,  as  before,  are  radial  elements  from  shroud  to  hub. 

The  thickness  of  the  blades  is  0.020  in.  at  the  inlet  and  is  variable  along 
*:he  passage  similar  to  impulse  blades,  thus  attaining  even  channel  cross- 
sectional  area  at  increased  curvature  section.  The  blades  are  tapered  toward 
the  hub. 


Diffuser  Design 


Designs  of  the  two  diffusers  were  completed  which  matched  the  impellers 
with  the  60°  backward  curved  and  30°  forward  curved  vanes.  A  vaned  diffuser 
was  selected  because  of  better  efficiency  and  smaller  overall  diameter  when 
compared  to  corresponding  vaneless  designs. 
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Peak  diffuser  performance  is  vital  especially  in  connection  with  the 
forward  curved  impeller,  because  a  higher  kinetic  energy  must  be  converted 
into  pressure  in  this  lower  reaction  case.  To  assize  a  maximum  performance, 
the  following  factors  must  be  incorporated: 

a.  Minimum  friction  (short  length  and  rapid  deceleration). 

b.  Minimum  boundary  layer  thickness  (short  'ength  with  slow  deceleration). 

c.  Retarded  separation  (slow  deceleration). 

d.  Maximum  diffusion  (rapid  deceleration  and  long  length). 

These  factors  are  somewhat  contradictory.  Therefore,  determir at  ion  of  an 
opt imum  comb inat i on  is  necessary. 

The  design  chosen  is  composed  of  th-ee  portions  ( refer"  to  Figure  20  ): 

a-  Minimum  vaneless  annulus  around  the  impeller. 

b.  Scroll  design  up  to  the  throat  area. 

c.  A  straight  diffusing  channel. 

The  optimum  recovery  of  a  straight  diffusing  channel  was  investigated  in 
Reference  7.  Figure  24  shows  the  influence  of  the  boundary  layer  thickness 
in  the  throat  on  the  pressure  recovery  of  the  diffuser  channel.  This  published 
curve  was  chosen  to  calculate  the  overall  efficiency  of  the  diffuser. 

The  dirfuser  inlet  portion  incorporating  the  vaneless  annulus  and  scroll 
contours  was  calculated  together.  The  basis  of  the  calculation  is  to  provide 
a  natural  deceleration  of  the  flow  induced  by  the  radius  increase  and  the 
friction  on  the  walls.  This  theory  is  assumed  to  produce  an  effective  deceler¬ 
ation  rate  and  low  friction  without  the  danger  of  separation  due  to  a  pressure 
gradient  that  is  too  large. 

The  basic  momentum  equation  that  follows  is  in  accordance  with  the  above 
theory  (designations  are  shewn  in  Figure  23  ): 


d(C  r)  f  C  r(2h  -  2h  +  b) 
dM  u  _ _ u  o 

dm  dx  ( h  -  h  )  (b  +  b) 

o  o 

and  the  continuity  equation: 


a*  V 

dh  -  :bo  CM(  r,:.1 


h  )  tan  y  +  b  +  b 
o  oj 
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where  M 


momentum 


m  =>  mass 

=  peripheraJ  velocity 
r  =  radius 

x  =  distance  in  velocity  direction 

C  =  velocity 

f  =  friction  factor 

h  =  free  channel  height 

b  =  channel  width 

=  radial  velocity 
V  =  angle  of  diffuser  side  wall 

The  two  differential  equations  above  were  integrated  and  the  resulting 
equations  were  programmed  on  the  IBM  computer.  The  computer  output  determined 
the  geometry  up  to  the  throat,  the  velocity  and  pressure  distributions  and  the 
losses  up  to  the  individual  points- 


Using  these  results,  the  boundary  layer  buildup  was  calculated  up  to  the 
throat  and,  with  this  value,  the  optimum  performance  (recovery)  of  the  straight 
diffuser  portion  was  estimated  in  accordance  with  Reference  7  and  Figure  24. 


Figure  33. 


SECTION  4 


MECHANICAL  DESIGN 


GENERAL 

As  was  the  case  in  the  hydrodynamic  design  of  the  pump,  mechanical 
design  was  strongly  affected  by  the  high  fluid  inlet  temperature,  and 
subsequent  selection  of  the  operating  fluid.  The  final  mechanical  arrange¬ 
ment  of  the  four- stage  pump  is  shown  in  Figure  34.  This  design  was 
based  originally  on  use  of  400°F  fluid  temperature,  and  was  later  modified 
to  accommodate  the  600°F  inlet  temperature. 

Design  problems  vie  re  anticipated  in  the  area  of  seals  and  bearings. 

A  seal  test  rig  was  built,  and  a  simple  controlled  gap  floating  carbon 
bushing  was  developed  to  provide  interstage  sealing.  Metallic  ”0"  rings 
were  used  for  static  seals,  and  bellows  type  carbon  face  seals  were  used 
for  shaft  seals. 

Silverplated  steel  sleeve  bearings  were  used  to  support  the  shaft, 
and  after  some  problems  with  a  hydrodynamic  pocket  bearing,  a  tilting  pad 
bearing  design  was  incorporated  fo~  the  thrust  bearing. 

The  unit  was  designed  to  be  shaft  driven  by  either  a  dynamometer 
(gear  box  input)  or  an  air  turbine  motor  (direct  driven^. 

FLUID  SELECTION 

The  specific  requirements  of  the  Problem  Statement,  Section  1, 

No.  1,  outlines  as  a  design  goal  that  the  pump  shall  handle  MLO  8200 
fluid  or  equivalent  from  room  temperature  to  600°F. 

Since  published  data  on  MLO  8200  fluid  (manufactured  by  Oronite 
Chemical  Company!  does  not  cover  the  temperature  range  above  500°F,  and 
for  some  characteristics  above  400°F,  the  Oronite  Chemical  Company  was 
contacted,  first  to  obtain  any  additional  data  on  MLO  8200  fluid  which 
might  be  available,  and  secondly,  to  obtain  data  on  equivalent  fluids. 

The  manufacturer  could  not  recorrmend  the  use  of  MLO  8200  for 
temperatures  above  520°F  and  suggested  the  use  of  Oronite  70  to  meet 
this  problem.  Data  comparing  the  properties  of  Oronite  70  with  MLO  820C 
is  appended  at  the  end  of  this  report. 

One  of  the  major  d i sadvantages  of  the  MLO  8200  fluid  is  the  decay 
in  viscosity,  at  elevated  tempe ratures ,  ’he  lower  limit  being  the  viscosity 
of  the  base  stock.  To  show  the  range  of  viscosities  expected,  an  extra¬ 
polated  curve  of  viscosity  as  a  function  of  temperature  for  both  the  base 
stock  and  MLO  8200  fluid  is  shown  in  Figure  35. 
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Figure  '5.  Viscosity  of  M  LO  8200  Fluid 
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In  spite  of  the  reported  disadvantages  of  using  MLO  8200,  the 
unit  was  designed  for  operation  in  this  fluid.  One  major  consideration 
was  the  opportunity  to  show  the  mechanical  effects  of  pump  operation 
in  a  fluid  with  obvious  poor  lubricity  and  lubrication  characteristics. 

BEARING  SELECTION 

Ball  Bearings 

For  this  application,  major  considerations  favoring  the  use  of 
ball  bearings  were  their  ability  to  take  axial  thrust  loads  and  to 
maintain  rotor  alignment.  Tests  conducted  at  AJResearch  indicated  ball 
bearings  of  the  required  size  could  be  operated  successfully  at  45, 000  rpm 
while  submerged  in  oil.  These  tests  also  showed,  however,  that  the 
resultant  bearing  losses  might  be  excessively  high  for  this  application. 

Testing  accomplished  by  both  New  Departure  and  Fafnir  Bearing 
Company  indicate  a  very  poor  life  for  204  size  bearings  operating  under 
moderate  loads  and  at  speeds  of  approximately  one  million  DN  (diameter 
x  rpm). 

The  following  is  a  summary  of  test  results  available  from  New 
Departure: 

1.  One  group  of  M-50  bearing  running  at  350°F  under  a  load 
factor  of  0.9  and  flooded  with  Oronlte  8200  ran  15  hours. 

On  examination,  they  were  found  to  be  in  excellent  condi¬ 
tion. 

2.  In  another  group  running  at  450°F  under  similar  conditions, 
all  failed  after  33  hours  operation. 

3.  A  viscosity  study  showed  that,  between  350°F  and  450°F 
and  over  a  time  span  of  150  hours,  Oronite  8200  under¬ 
goes  a  severe  viscosity  reduction.  This  was  found  to  be 

*  as  high  as  40  percent.  Due  to  this  reduction,  the  lubricant 

is  capable  of  supporting  only  extremely  light  loads.  Mod¬ 
erate  loads  (producing  load  factors  of  1.0  and  3.0)  result 
in  metal  to  metal  contact  and  eventual  seizure.  Heavy  loads 
(producing  load  factors  up  to  0.9)  result  in  bearing  life  in 
the  range  of  I  percent  of  calculated  B-10  life. 

Based  on  the  above  factors  and  the  fact  that  there  Is  further 
viscosity  reduction  in  the  600°F  vicinity,  it  was  recommended  that,  to 
obtain  any  degree  of  bearing  life,  the  loads  must  be  extremely  light 
(load  factors  above  5.0).  For  the  204  size,  the  size  contemplated  for 
this  application,  operating  at  50,000  rpm  and  at  600°F  to  650°F  under 
moderate  loads,  it  would  be  impractical,  if  not  impossible,  to  project 
any  life  forecast.  New  Departure  engineers  predicted  it  to  be  in  the  area 
of  10  hours  or  less. 
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Because  of  this  evaluation  and  other  similar  information,  it  has 
been  decided  to  use  journal  bearing  configurations  exclusively. 

Sleeve  Bearings 

Suitable  materials  for  sleeve  bearings  such  as  graphite, 

Meehanite  AQ  and  silver  or  silver  plated  steel  are  available.  Sliver 
plated  steel  journals  were  chosen.  These  bearings  are  compatible  with 
the  fluid,  have  sufficient  imbeddibi 1 i ty  to  accommodate  some  fluid 
contamination,  and  are  relatively  simple  to  fabricate. 

Because  inlet  pressure  of  100  psi  was  available  as  a  design  cri¬ 
terion,  use  of  a  simple  hydrostatic  bearing  was  selected.  The  inboard 
journal  was  lubricated  by  moderately  high  pressure  fluid  bled  from  the 
first  stage  discharge,  and  the  outboard  journal  was  lubricated  by  fourth 
stage  discharge  fluid  leakage. 

Thrust  Bearings 

Because  sleeve  bearings  were  selected  for  the  final  design,  a 
thrust  bearing  was  required  to  locate  the  rotor  axially  and  absorb 
any  unbalanced  hydraulic  thrust  load.  Five  thrust  bearing  designs 
were  investigated  using  a  preliminary  figure  for  unbalanced  thrust 
load.  Of  these,  three  appeared  to  be  satisfactory.  The  first  design 
is  a  tilting  pad  thrust  bearing  manufactured  by  Industrial  Tectonics, 

Inc.  This  is  a  self-leveling  unit  and  is  unique  in  that  the  rotating 
thrust  collar  is  graphite  and  the  pads  are  steel  (the  reverse  of  the 
usual  construction).  Since  graphite  is  compatible  with  MLO  8200  fluid 
and  is  suitable  for  use  at  the  temperatures  involved,  it  was  satisfactory 
for  the  appl icat ion. 

The  second  design  considered  acceptable  is  the  stepped  pad  con¬ 
figuration  which  has  the  same  normal  load  carrying  capacity  as  a  tilting 
pad  bearing  and  approximately  the  S' me  losses.  It  is  not  self-aligning 
and  probably  does  not  have  quite  as  much  reserve  load  carrying  ability 
for  a  given  design  load  as  the  tilting  pad  bearinq.  However,  it  is 
simpler  mechanically  and  relatively  easy  to  manufacture. 

The  third  design  considered  applicable  is  the  hyd>*odynam: c  pocke^ 
thrust  bearing.  This  unit  is  more  difficult  to  manufacture  than  the 
stepped  pad  bearing  because  of  close  tolerances  even  though  it  is  simple 
mechanically.  Its  losses,  on  the  same  design  basis,  are  approximately 
50  percent  of  the  losses  of  either  the  tilting  pad  bearing  or  the 
stepped  pad  bearing. 

The  initial  four  stage  pump  incorporated  the  hydrodynamic  pocket 
bearing  design.  This  design  was  chosen  on  the  basis  of  the  extremely 
low  bearing  losses.  However,  initial  runs  using  MIL  L-6081  fluid 
indicated  that  extensive  development  would  be  required  in  order  to  fully 
utilize  the  advantages  of  the  pocket  bearing.  Therefore,  the  design  was 
modified  to  use  the  tilting  pad  design. 
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The  Industrial  Techtonics  standard  bearing  was  modified  to  use  a 
thrust  runner  Hide  from  hardened  17-4  PH  steel.  The  thrust  surfaces 
were  nitrided  to  provide  an  extremely  hard  wearing  surface.  The 
standard  steel  shoes  (tilting  pads)  were  silver  plated  in  the  same 
manner  as  the  journal  bearings. 

The  tilting  pad  thrust  bearing  develops  its  own  hydrodynamic 
pressure,  and  requires  lubrication  only  for  cooling. 

LUBRICATION 

In  order  to  provide  a  self-contained. uni t,  bleeds  were  provided 
at  each  of  the  intermediate  stages  to  provide  bearing  lubrication  to 
both  sides  of  the  thrust  bearing  and  to  the  inboard  journal  bearing. 
Throughout  the  testing,  only  first  stage  bleed  was  used.  The  bleed 
pressure  *as  limited  to  approximately  200  psi  maximum  and  bearing  flows 
were  kept  at  approximav -»lv  0.5  gallon  per  minute  per  bearing. 

The  outboard  journal  bearing  was  lubricated  by  leakage  past  the 
fourth  stage  seal.  During  preliminary  room  temperature  runs,  this 
leakage  flow  was  measured,  and  results  were  in  good  agreement  with 
values  expected  from  seal  leakage  curves. 

SEAL  SELECTION 

The  following  static  seals  were  considered  for  this  application: 
elastomer  "0"  rings,  metallic  "0"  rings,  crush  gaskets,  and  proprietary 
metallic  face  seais. 

Elastomer  "0"  Rings 

A  careful  study  of  the  various  elastomer  materials  presently  avail¬ 
able  has  revealed  that  there  are  none  now  available  which  are  both 
compatiMe  with  ML0  8200  fluid  and,  at  the  same  time,  suitable  for 
operation  at  600°F  for  periods  in  excess  of  a  few  hours.  Consequently, 
this  approach  was  abandoned. 

Metallic  "0"  Rings 

A  review  of  literature  and  manufacturer  catalogs  indicated  a  vented 
silver-plated  stainless  steel  "0"  ring  ciuld  be  successfully  utilized  in 
this  application.  These  seals  depend  on  flange  loading  to  deform  the  seal, 
however,  in  this  development  unit,  heavy  flange  sections  presented  no 
problems. 

Crush  Gaskets  *-'*  , 

■  ''  '■■ii  ■■  . 

* 

Asbestos-filled  gaskets,  such  as  "Flex i tal 1 ic",  manufactured  by 
Flexitallic  Gasket  Company  and  "Guardian",  manufactured  by  Garlock,  Inc., 
are  spirally  wound  gaskets  which  have  been  used  successfully  for  years  for 
high  pressures  at  temperatures  considerably  higher  than  required  by  this 
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application.  They  do  not  require  surfaces  finished  as  highly  as  other 
seals,  and  are  relatively  inexpensive.  However,  they  require  high 
compression  loading  to  be  effective.  There  is  no  published  data  on 
loading  requirements,  nor  any  basis  for  design  studies,  since  these 
gaskets  are  supplied  for  standard  flanges.  However,  discussions  with 
one  supplier  indicate  the  compressive  loading  should  be  taken  as 
approximately  4.7  times  the  pressure  to  be  sealed.  Even  with  a  very 
narrow  gasket,  this  required  a  much  higher  loading  than  the  published 
figures  for  vented  metallic  "0"  rings  for  the  same  application.  In 
view  of  this,  investigation  of  this  type  of  sea!  was  abandoned. 

Proprietary  Metallic  Face  Seals 

4  number  of  proprietary  metallic  face  seals  were  investigated  for 
this  application.  From  testing  done  for  other  applications,  there  was 
no  question  that  several  of  these  were  satisfactory  for  use  in  the  pump. 
They  all  require  exceedingly  well  finished  contact  surfaces  to  be 
effective.  These  surfaces  can  be  easily  damaged  during  assembly  and 
disassembly.  Since  this  was  a  development  unit,  frequent  assembly  and 
disassembly  was  anticipated  during  the  course  of  testing.  Consequently, 
it  was  felt  that  these  seals  should  not  be  used  unless  unexpected 
problems  developed  with  metallic  "0”  rings. 

Dynamic  Shaft  Seals 


Since  the  four-stage  pump  configuration  would  have  the  first- 
stage  inlet  at  the  end  where  the  shaft  emerges,  the  highest  pressure  to 
be  sealed  against  atmosphere  on  the  rotating  shaft  was  a  100  psig 
maximum  inlet  pressure.  A  face  contact  sea!  with  a  metallic  bellows 
secondary  seal  was  used.  The  pressures,  temperatures,  and  rubbing  speeds 
involved  were  well  within  the  present  state  of  the  art.  One  question 
raised  was  the  exact  state  of  the  seepage  residue  of  600°F,  MLO  8200 
fluid  after  exposure  to  the  atmosphere.  However,  the  problem  was  not 
considered  to  be  severe  enough  to  prevent  the  use  of  these  seals. 

Impeller  and  Interstage  Shaft  Seal 

The  high  pressure  rise  per  stage  of  this  pump,  the  high  rotative 
speed  and  the  small  space  available  all  combine  to  render  the  impeller 
and  inter-stage  shaft  seals  a  very  difficult  problem.  Several  approaches 
were  considered.  Among  them  were  face  contact  seals,  labyrinth  seals, 
and  control led-gap  floating  bushing  seals. 

1 .  Face  Contact  Seals 


This  type  of  seal  has  desirable  characteristics  in  that  leakage  and 
friction  losses  should  be  minimal;  however,  the  problems  of  application 
to  this  pump  are  quite  formidable..  The  lack  of  a  suitable  elastomer  for 
600°F  service  for  use  as  a  secondary  seal  limited  the  selection  to 
metallic  bellows  type  secondardy  seals.  Several  manufacturers  of  face 
contact  seals  were  contacted  and  two,  Koppers  Co.  and  Sealol,  offered 
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a  double  bellows  supported  face  contact  seal.  The  space  required  for 
this  arrangement  was  far  in  excess  of  that  available.  This  factor,  plus 
the  problem  of  assembly  and  maintenance  of  the  required  normality  and 
flatness  of  the  mating  rotating  surface,  caused  abandonment  of  this 
approach. 

2.  Labyrinth  Seals 

Preliminary  calculations  indicate  that  a  radial  clearance  sub¬ 
stantially  less  than  0.001  inch  must  be  maintained  to  reduce  leakage 
losses  of  a  labyrinth  seal  to  an  acceptable  value.  These  close 
clearances  might  be  maintained  in  service  using  a  ball  bearing 
supported  rotor.  However,  since  sleeve  bearings  were  selected  for  the 
final  design,  labyrinth  seals  could  not  hold  close  radial  clearances 
requi red. 

3.  Control led-Gap  Floating  Bushing  Seals 

This  type  of  sea!  limits  the  flew  between  the  rotating  and  stationary 
elements  by  maintaining  a  very  close  radial  clearance  between  the 
rotating  element  and  a  non-rotating  floating  bushing.  The  secondary 
seal  is  obtained  by  the  pressure  unbalance  holding  the  bushing  axially 
against  a  radial  lapped  surface  on  the  housing. 

Three  different  designs  were  considered.  All  three  can  be  accom¬ 
modated  within  the  available  space,  and  should  have  approximately 
equivalent  leakage  and  friction  losses.  They  differ  only  in  materials 
and  method  of  maintaining  the  desired  clearance  over  the  operating 
temperature  range. 

The  first,  proposed  by  Chicago  Rawlvde  Manufacturing  Co.,  consists 
of  a  graphite  bushing  shrunk  in  a  bushing  of  the  same  materia!  as  the 
rotating  element.  The  graphite  is  compatible  with  MLO  8200  fluid  and 
Is  also  a  good  bearing  material  at  the  highest  anticipated  temperature. 

The  outer  bushing  serves  to  control  the  thermal  expansion  rate  of  the 
inner  carbon  bushing  which  is  always  under  compression.  Since  the 
outer  ring  is  of  the  same  material  as  the  rotating  element,  the  thermal 
expansion  rates  would  be  the  same  and  the  clearance  between  rotating 
element  and  stationary  floating  bushing  would  remain  constant  over  the 
range  of  temperatures  expected. 

The  second  arrangement,  proposed  by  Koppers  Company,  consists  of 
a  stationary  floating  bushing  of  440  C  steel  hardened  to  a  60-57  Rockwell 
''C"  scale  hardness  mated  against  a  hard  chromeplated  surface  on  a 
17-4PH  steel  rotating  element.  This  arrangement  would  show  a  slight  decrease 
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in  radial  clearance  with  increasing  temperature,  a  desirable  feature  when 
properly  controlled. 

The  third  arrangement  consisted  of  a  stationary  floating  bushing  of 
Meehanite  AQ  which,  in  the  temperature  range  under  consideration,  would 
have  the  same  temperature  coefficient  of  expansion  as  1 7—4 PH  steel,  the 
material  used  for  the  rotating  element. 

A  seal  test  rig  suitable  for  evaluating  control  I ed-gap  float, ng  bushing 
seals  was  designed  and  fabricated.  A  sketch  of  this  test  rig  is  included  in 
Figure  36.  Figures  37  and  38  are  photographs  of  the  complete  test  setup 
including  the  air  turbine  drive  motor.  Figure  39  depicts  the  Meehanite  AQ 
version  of  an  impeller  and  interstage  shaft  seal.  Alternate  versions  of 
this  seal  include:  (l)  A  I7-4PH  stainless  steel  floating  bushing  with  a 
silver  plated  bore,  (2)  A  carbon  bushing  shrunk  in  a  I7-4PH  ring  (Chicago 
Rawhide  Mfg.  Co.),  (3)  A  440C  stainless  steel  ring  mated  with  a  I7-4PH  chrome 
plated  shaft  (Koppers  Co.,  Inc.),  and  (4)  A  carbon  bushing  shrunk  in  an 
AM-355  stainless  steel  ring  mated  with  a  I7-4PH  chrome  plated  shaft  (Koppers 
Co. ,  Inc. ). 

The  seal  test  rig  was  arranged  to  measure  rotational  speed,  the  leakage 
rate  past  the  seals  and  power  input.  The  air  turbine  drive  motor  was 
selected  for  this  application  because  of  its  availability  and  flexibility 
of  speed  and  power  control.  The  test  setup  had  been  arranged  such  that 
tests  were  conducted  using  a  variety  of  fluids  including  low  viscosity  oil 
which  is  equivalent  in  viscosity  at  moderate  temperatures  to  MLO  8200  Fluid 
at  high  temperatures. 

Results  of  this  test  program  are  presented  in  Section  5  •  As  a  direct 
result  of  these  tests,  interstage  seal,  similar  to  item  (2)  above  was 
successfully  developed  and  incorporated  into  the  design  of  the  pump.  One 
change  in  the  seal  design  made  as  a  result  of  the  experience  gained  on  the 
single-stage  test  rig  was  the  addition  of  a  wave  spring  to  provide  a  posi¬ 
tive  axial  seating  force  at  the  seal.  This  design  change  eliminated  the 
tendency  of  the  seal  to  not  seat  as  the  pump  was  started. 

BEARING  CRITICAL  SPEEO  CONSIDERAl IONS 

A  dynamic  analysis  of  the  bearing  system  was  made  to  investigate  possible 
bearing  critical  speed  problems.  This  analysis  indicated  the  existence  of 
the  possibility  that  operation  at  45,000  rpm  might  be  close  to  a  calculated 
critical  speed.  The  analysis  was  dependent  on  assumptions  made  regarding 
relative  shaft  diameters  and  the  spring  stiffness  of  the  bearings  at  dif¬ 
ferent  speeds  and  under  varying  load  conditions.  Results  of  this  analysis 
are  presented  in  Figure  ,40,  which  shows  bearing  spring  rate  versus  shaft 
speed  at  a  600°F  fluid  temperature.  This  figure  also  depicts  critical  speed 
modes  at  these  spring  rates. 

It  was  felt  that  even  if  a  critical  speed  coincides  with  operating  sp-ed, 
the  hydraulic  damping  afforded  by  having  the  complete  rotating  assembly  sub¬ 
merged  in  the  hydraulic  fluid  would  prevent  excessive  shaft  excursion 
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Figure  36.  Dynamic  Seal  Test  Rig  (Cross  Section) 


Figure  37.  Seal  Test  Setup 
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38.  Seal  Test  Setup 
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Figure  AO.  Cnticai  Speed  Analysis,  A  SO  Hydrauiic  P^o 
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Actual  testing  of  the  four  stage  pump  indicated  that  this  damping  did 
occur  as  no  shaft  criticals  or  excursions  were  noted  at  any  ooerating  speed- 
up  to  47,000  rpm. 

ASSEMBLY  f^OCEDU  c 

All  f>  ur  pump  impellers  were  stacked  on  che  drive  shaft  and  locked  axially 
a~ainst  a  shoulder  by  means  of  a  locxing  nut.  Positioning  of  the  impellers 
with  respect  to  their  mating  diffusers  and  maintenance  of  proper  shoulder 
clearance  was  accomplished  by  a  combination  of  shimming  of  the  diffusers, 
and  gi Inding  the  impeller  hubs.  As  subsequent  ass  mlies  and  disassemblies 
were  made,  it  became  increasingly  difficult  to  maintain  alignment.  This 
fact  became  apparent  during  actual  testing  because  of  the  low  output  flows 
attributed  to  this  problem. 

Figures  4*  and  42  show  the  detail  part'  of  the  pu~o.  Figure  43 
shows  a  parts  a.  ray  of  a  disassembled  pump  ano  Figure  44  shows  the  com¬ 
plete  assembled  unit. 

TURBINE  DRIVE 

In  order  to  achieve  the  smoothest  possible  driving  mechanism  for  the 
initial  break-in,  i  <.  was  decided  to  adaot  an  ex  sting  air  turbine  drive 
to  the  pump.  The  compressor  section  ot  cabin  air  compressor  was 

discarded  and  a  special  adapter  was  mad  o  match  the  turbine  housing  and 
pump  inlet  flange  bolt  circles.  The  soiii  »  was  the  same  as  that  used  on 
the  pump  so  that  no  additional  rework  was  needed. 

Speed  control  was  accomplished  by  va  /ing  turbine  inlet  air  pressure, 
and  by  varying  turbine  nozzle  area.  standard  pneumatic  nozzle  positioner 
was  used  for  this  latter  purpose  in  conjunction  with  pressure  regulated  a;  r 
supply. 

Using  a  turbine  to  p  Dvioe  a  direct  driving  means  for  the  pump  eliminates 
the  need  for  a  gear  box  and  eliminates  the  shock  of  step  speed  charges. 
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SECTION  5 


EXPERIMENTAL  TESTING 


GENERAL 

This  section  is  divided  into  three  parts.  The  first  is  devoted  to 
individual  component  testing  of  seals  and  a  single  stage  pump,  the  second 
to  the  final  pump  design  and  the  i ns t rumen  La  lion  provided  and  the  third 
part  to  the  discussion  of  pump  testing  and  results  obtained. 

COMPONENT  TESTING 

Prior  to  the  finalization  of  the  design  of  the  four  stage  pump, 
extensive  development  testing  was  accomplished,  primarily  of  the  inter- 
stage  seal,  and  the  hydrodynamic  performance  of  the  calculated  impeller 
and  diffuser  design. 

These  component  tests  resulted  in  the  development  of  an  interstage 
seal  suitable  for  this  application  and  sufficient  data  was  obtained  to 
permit  engineering  application  of  this  seal  in  other  designs. 

Design  of  the  hydrodynamic  components,  as  predicted  by  the  computer 
program  was  confirmed  by  means  of  testing  accompli  hed  on  a  single  stage 
rig.  In  addition  to  design  confirmation,  these  tests  permitted  further 
use  of  the  computer  program  with  an  increased  confidence  factor. 

Interstage  Seal  Tests 

Testing  of  various  seal  configurations  resulted  in  the  selection  and 
development  of  a  controlled  gap  floating  carbon  bushing  seal  manufactured 
by  Chicago  Rawhide  Company,  Part  Number  800561.  This  seal  is  shown  in 
Figures  45  and  46.  Testing  was  done  on  the  air  turbine  seal  test  rig 
shown  in  Figures  47  and  48.  Test  runs  were  conducted  at  various  speeds, 
fluid  viscosities  and  pressure  differentials. 

Initial  testing  using  the  seal  test  rig  resulted  in  two  cases  of  shaft 
excursion  which  caused  seal  damage.  This  was  investigated  by  the  use  of 
inductive-type  proximity  pick-ups  and  found  to  be  caused  by  a  shaft  critical 
speed  at  about  41,000  rpm.  The  problem  was  solved  by  reducing  the  weight  of 
the  rotating  parts  on  tkD  seal  end  of  the  unit  which  increased  the  critical 
speed  to  above  the  45,000  rpm  operating  speed. 

A  thorough  investigation  of  the  Chicato  Rawhide  Seal  800561  operating 
with  a  radial  clearance  of  0.0005  in.  was  conducted.  Test  results  indicated 
the  leakage  at  the  viscosity  corresponding  to  ML0  8200  at  600°F  (1.98 
centistokes)  was  0.29  gpm  at  a  sea!  pressure  differential  of  750  ps i  which  is 
very  c'ose  to  the  leakage  assumed  in  the  original  pump  design  calculations. 
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Figure  45.  Interstage  Seal 
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This  initial  testing  was  conducted  using  jet  engine  lubricating  oil 
(H1L-L-6O0IB,  Grade  1005).  Leakage  values  were  measured  at  various 
constant  speeds  as  shown  on  Figure  49  which  indicate  the  effect  of 
speed  on  the  leakage  was  slight.  Leakage  at  the  45,000  rpm  speed  was 
measured  at  three  pressure  differentials,  750  psi,  100  psi  and  1250 
psl  as  shown  on  Figure  50.  The  seals  are  expected  to  operate  at 
200  to  1000  psi  differentials  in  the  pump.  In  order  to  invest! gage  the 
operation  of  the  seals  at  higher  fluid  viscosities,  a  change  to  Sinclair 
L-1048  was  made  and  data  obtained  at  viscosities  up  to  29.5  cent (stokes 
as  shorn  In  Figure  51.  Operation  of  the  seals  up  to  this  point  was 
normal,  although  the  oil  temperature  rise  became  quite  large  (the 
temperature  differential  through  seal  at  the  29.5  ccntistokes  condition 
was  I73°F).  At  the  56  centistokes  condition  after  about  three  minutes 
run  time,  the  leakage  rate  suddenly  increased  and  the  oil  temperature 
rise  decreased.  On  disassembly  of  the  unit,  it  was  found  that  the  seat 
had  rubbed  on  the  downstream  half  of  the  bore  only  which  resulted  in 
increased  clearance  and  consequently,  the  increased  leakage. 

It  appears  that  this  rubbing  on  one  portion  of  the  seal  was  caused 
by  expansion  of  the  rotating  sleeve  in  this  local  area  due  to  the  consid¬ 
erable  oil  temperature  rise  as  It  passed  through  the  seal.  In  the  config¬ 
uration  tested  in  the  seal  test  rig,  there  appears  to  be  a  viscosity  limit 
of  approximately  50  centistokes  maximum  with  a  radial  clearance  of  0.0005 
inch  at  45,000  rpm  and  1000  psi  pressure  differential.  However,  it  was 
felt  that  in  the  actual  pump  configuration,  the  viscosity  limit  would  be 
considerably  higher  due  to  better  heat  transfer  from  the  seal  areas  on 
the  impeller  to  the  large  volume  of  fluid  passing  through  the  pump.  For 
comparison,  the  viscosity  ML0  8200  hydraulic  fluid  at  70°F  is  49  centi¬ 
stokes  and  at  I07°F  is  30  centistokes. 

Measurements  were  made  of  the  power  losses  occurring  in  the  seals, 
by  measuring  the  power  supplied  by  the  air  turbine  (using  turbine  air¬ 
flow  and  turbine  air  temperature  drop  data  and  deducting  turbine  bearing 
power).  Also,  power  losses  were  determined  by  measurements  of  the  leakage 
flow  rate  and  temperature  rise  of  the  leakage  oil  as  it  passed  through  the 
seal.  Plots  of  power  losses  based  on  the  latter  data  are  shown  in  Figures  52, 
53  and  54.  The  power  values  measured  by  the  two  methods  agreed  quite 
well  at  viscosity  values  up  to  8  centistokes.  Above  this  point,  the  power 
measured  b*/  turbine  air  horsepower  became  higher  than  the  oil  horsepower 
values  due  to  the  increased  drag  losses  of  the  rotating  parts  at  locations 
other  than  the  seals  operating  in  high  viscosity  oil. 

Throughout  the  calibration  of  the  60  degree  backward  curved  impeller 
in  the  single-stage  test  rig,  the  floating  bushing  seals  and  carbon  face 
type  shaft  seal  performed  very  well.  The  leakage  of  the  outboard  bushing 
seal  was  measured  throughout  the  test.  The  leakage  rate  agreed  closely 
with  the  measurements  made  on  the  seal  test  rig,  ranging  from  0.143  gpm 
to  0.261  gpm  at  the  45,000  rpm  speed.  The  internal  pressure  measuring 
ports,  which  were  incorporated  in  the  test  rig,  permitted  monitoring  of 
the  impeller  inlet  eye  seal  pressure  differential  although  the  leakage 
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rate  could  not  actually  be  measured  at  this  location.  The  pressures  on 
each  side  of  the  impeller  were  nearly  the  same.  This  indicates  similar 
leakages  through  both  seals  and  aiso  means  that  the  axial  thrust  was 
close  to  being  balanced.  The  leakage  measurements  and  internal  pressure 
measurements  did  not  indicate  any  increase  in  seal  clearances  during  the 
test  run. 

Single  Stage  Impeller  Tests  ' 

To  confirm  the  results  uf  the  analytical  hydrodynamic  design,  a 
singi-*  stage  impeller,  incorporating  60  degree  backward  curved  blades 
was  fabricated  and  tested.  Figures  55  and  56  show  this  impeller.  Test  rig 
details  are  shown  in  Figures  57  through  61,  and  the  test  setup  is 
shown  in  Figures  62  through  66  .  The  test  rig  was  driven  by  means 

of  a  relatively  low  speed  dynamometer  and  high  speed  gear  box.  This 
gear  box  is  identical  to  that  used  on  AiResearch  Compressor  Part  206420 
(Lockheed  L- 1 88  Cabin  Compressor). 

All  testing  on  the  single  stage  rig  was  accomplished  using 
MIL-L-608IB,  Grade  1005  lubricating  oil  at  I 80°F  to  I90°F.  This  results 
in  a  kinematic  viscosity  of  approximately  2  centistokes,  the  same  as  that 
of  ML0  8200  fluid  at  600°F. 

Results  of  these  tests  are  shown  in  Figures  67  and  68  .  Figure  67 
shows  pressure  rise  versus  flow  at  constant  speed  parameters  from  15,000  to 
45,000  rpm.  Constant  efficiency  parameters  are  super- imposed  showing  a  peak 
efficiency  island  of  65  percent  of  approximately  40,000  rpm  and  efficiency 
of  64.5  percent,  mcximum,  at  the  maximum  speed  of  **5,000  rpm.  Shaft  power 
input  versus  flow  is  presented  in  Figure  for  the  same  constant  speed 

parameters. 

Data  shown  at  a  speed  of  42,000  rpm  is  particularly  significant 
because,  when  operating  at  this  speed,  the  pump  pressure  rise  is  exactly 
equal  to  the  pressure  required  in  the  four  stage  4000  psi  pump  at  the 
design  point  (using  MIC  8200  fluid  at  600°F).  The  efficiency  shown  on 
these  curves  is  an  overall  efficiency  based  on  hydraulic  power  output  and 
pump  shaft  power  input.  Thus,  it  includes  power  losses  in  two  ball  bear¬ 
ings  (runnirg  in  an  air-oil  mist),  one  carbon  face  shaft  seal  and  two 
floating  bushing  seals  (one  at  the  impeller  inlet  eye  and  one  on  the  out¬ 
board  side  of  the  wheel).  The  shaft  nower  input  was  measured  by  an  electric 
motor  dynamometer  with  deduction  of  power  losses  in  the  12.26:1  ratio  step-up 
qear  box.  These  gear  box  losses  were  based  on  measurement  of  the  gear  box 
lubricating  oil  flow  rate  and  oil  temperature  rise  through  the  gear  box. 

An  investigation  was  conducted  to  determine  the  effect  of  increased 
fluid  viscosity  on  pump  overall  efficiency.  Performance  was  measured  at 
the  42,000  rpm  speed  at  viscosities  of  2,  3  and  4  centistokes.  This  was 
accomplished  by  cooiing  the  oil  to  increase  its  viscosity.  It  was  found 
that  there  was  a  reduction  in  peak  efficiency  of  approximately  1.5  percent 
for  each  centistoke  increase  in  viscosity.  These  results  are  shown  in 
F i qu  re  69  . 
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Figure  55.  Backward  Curved  I  ">eMer 
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Figure  56.  60°  Backward  Curved  Impeller 
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57.  Single  Stage  Impeller  Test  Rig 


Figure  64.  Single  Stage  Test  Setup 
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Figure  68.  Input  Power  vs  Flow,  Single-Stage  Centrifugal  Pu»,p  Assembl 
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Single-Stage  Centrifugal  Pump  Assemb) 


The  performance  of  the  single-stage  backward  curved  impeller  was  quit* 
encouraging.  Calculations  made  using  this  test  data  show  an  Impeller  hydraulic 
efficiency  of  80  percent.  This  is  an  excellent  efficiency  for  an  impeller  of 
this  size  and  is  5  to  8  per  •’nt  higher  than  the  predicted  impeller  efficiency 
used  in  the  original  design  work. 

As  can  be  seen  from  the  pressure  rise  versus  flow  curves,  the  pump  exceeds 
the  required  pressure  rise.  If  the  pressure  rise  of  approximately  1300  psi, 
which  occurs  at  45,000  rpm,  is  corrected  for  the  difference  in  densities  cf 
ML0  8200  fluid  at  600°F  and  MIL-L-6081,  Grade  1005  at  I85°F  (MIL-L-6081  is 
14  percent  more  dense),  it  becomes  1148  psi  pressure  rise  for  ML0  8200  at 
600°F.  Since  the  required  pressure  rise  per  stage  is  975  psi,  the  impeller 
diameter  could  be  reduced  slightly  which  might  reduce  internal  disk  friction 
losses  thus  increasing  the  overall  efficiency  to  66  to  70  percent. 

The  next  test  conducted,  therefore,  was  a  performance  check  at  45,000  rpm 
with  2.3  inches  outside  diameter  backward  curved  impeller  trimmed  to  2.160  inches 
outside  diameter.  The  2.160  inch  diameter  was  selected  to  provide  the  pressure 
rise  required  for  the  four-stage  pump  at  45,000  rpm.  For  this  test,  no  changes 
were  made  in  either  the  diffuser  or  in  the  shrouds  in  order  to  expedite  perform¬ 
ance  of  the  test.  It  was  found  that  the  pump  efficiency  decreased  about  2  per¬ 
cent.  The  decrease  in  pump  pressure  rise  was  more  than  would  be  predicted  using 
an  !mpel ler  diameter  squared  relationship.  Pressure  rise  dropped  from  1300  psi 
to  1 i 00  ps i  and  the  stable  flow  range  of  the  pump  decreased  by  4  gpm.  However, 
it  is  probable  that  changing  the  diffuser  and  shroud  designs  to  match  the 
reduced  impeller  diameter  would  improve  performance.  The  final  four-stage 
pump  will  incorporate  the  2.3- inch  diameter  impeller. 

In  addition  to  the  60  degree  backward  curved  blade  impeller,  it  was 
originally  intended  to  conduct  a  series  of  similar  tests  incorporating  a 
30  degree  forward  curved  impeller.  The  impeller,  as  shown  in  Figure  69 
was  fabricated  and  initial  test  runs  started.  It  became  apparent  that  use 
of  this  impeller  caused  static  pressure  variations  between  the  impeller 
shroud  and  housing  which  lead  to  high  axial  thrust  loads  on  the  seals  and 
bearings.  Rather  than  develop  a  test  rig  to  overcome  this  problem,  it  was 
decided  to  concentrate  efforts  on  the  complete  four  stage  pump  which  was 
ready  at  the  same  time. 

FOUR- STAGE  PUMP 

Test  Setup 

Initial  test  break  in  runs  were  made  using  lubricating  fluid  per 
MIL-L-608 I B,  Grade  1005.  These  tests  w*re  conducted  at  ambient  room 
temperatures  with  a  fluid  irlet  temperature  of  approximately  I90°F. 

Figures  71  and  72  show  the  pump  test  setup  for  this  run. 
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Fiqure  70.  30°  Forward  Curve  Impel 
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Figure  71 


View  of  P uffp  and  TurDine  Drive  Test  Setup  tor  Run 
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Figure  72.  Pump  ana  Turbine  Drive  Test  Setup  for  Run  I 
(Viewed  from  Pump  End) 
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In  order  to  achieve  the  smoothest  possible  driving  mechanism  for  the 
initial  break  in,  it  was  decided  to  adapt  an  existing  air  turbine  drive 
to  the  pump.  The  compressor  section  of  an  AiResearch  Part  206170  (Boeing 
707)  Cabin  Air  Compressor  was  removed  and  a  special  adapter  was  made  to 
match  the  turbine  housing  and  pump  inlet  flange  bolt  circles-  Speed  control 
was  accomplished  by  varying  turbine  inlet  air  pressure,  and  by  varying  turbine 
nozzle  area.  The  standard  pneumatic  nozzle  positioner  was  used  for  this  purpose 
and  a  pressure  regulated  air  supply.  Figures  73  and  74  show  the  pump, 
adapter,  and  turbine  drive. 

Use  of  an  air  turbine  to  provide  a  direct  driving  means  for  the  pump 
eliminates  the  need  for  a  gear  box,  and  also  eliminates  the  shock  of  step 
speed  changes.  Measurement  of  pump  inlet  power,  however,  depends  on  an 
accurate  measure  of  turbine  air  weight  flow  and  temperature  drop.  Bearing 
losses  for  the  turbine  drive  have  been  well  established  by  many  previous 
tests . 
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Figure  75  presents  a  schemat’-  o,  he  setup  for  the  high  temperature 
tests  of  the  four  stage  hydrai  l'  pump  usin<“  MLO  3200  fluid.  Figures  76 
and  77  show  the  actual  test  installation. 

The  pump  generated  neat  raised  tie  fluid  temperature  to  the  des i red 
level,  and  this  temperature  was  maintained  by  bypassing  a  portion  of  the 
discharge  flow  through  a  heat  exchjnaer.  Bearing  leakage  flow  was  returned 
to  the  inlet  rese  voir  through  anothc~  heat  exchanger.  The  fluid  reservoir 
was  remotely  placed  in  nrde  to  maintain  the  bulk  temperature  of  the  fluid 
at  a  reasonably  iew  /a)ue  ( approx imatel y  400°F).  Lubricating  fluid  for  the 
thrust  bearing  and  the  irbcarc  journal  bearing  was  bled  from  the  first  stage 
di scharge. 

Because  of  safetv  cons 'derations,  the  pump  and  heat  exchanger  were 
blanketed  in  a  C02  atmosphere. 

Ins  trumentat i on 

During  break  in  runs  at  low  fluid  temperatures  (approximate^1  !90°F) 
and.  initial  high  temperature  runs,  shaft  deflection  instrumentation  was 
included.  Although  this  type  instrumentation  will  not  provide  valid  results 
at  elevated  temperatures,  sufficient  information  was  gathered  to  confirm  that 
no  excessive  shaft  excursions  had  occurred  during  the  runs. 

Total  lubricating  flow  to  the  thrust  bearing  and  inboard  journal  bearing 
was  measured  by  means  of  a  high  temperature  flow  sensing  element.  This  was 
used  in  conjunction  with  a  remotely  operated  needle  valve  to  regulate  and 
monitor  the  fluid  quantity  bled  from  the  first  stage. 

Outboard  journal  bearing  flow  vas  measured  during  the  room  temperature 
runs.  This  flow  was  in  close  agreement  v/ith  the  leakage  rates  as  measured 
or,  the  seal  test  rig,  indicating  a  bearing  flow  of  approximately  1.28  gallons 
per  minute  at  4000  ps i  discharge  pressure.  During  high  temperature  runs,  this 
flow  was  retuined  directly  to  the  inlet  reservoir  in  a  closed  loop  and  was  not 
measured. 
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Figure  73.  Four-Stage  Pump  Assembly  and  Turbine  Drive 
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Figure  74.  Four-Stage  Pump  Assembly  and  Turbine  Drive 
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Figure  76.  Test  Installation  for  the  High  Temperature  Tests  of  the 
Four-Stage  Hydraulic  Pump  (L208043)  Using  MIL-0-8200 
Hyd rau lie  Fluid 
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Figure  77.  Test  Installation  for  the  High  Temperature  Tests  of  the  Four-Stage 
Hydraulic  Pump  (L208043)  Using  MIL-0-8200  Hydraulic  Fluid 


Pump  oil  flow  was  measured  by  a  high  temperature  flow  sensing  element 
in  the  pump  inlet.  Oil  inlet  and  discharge  temperature  was  measured  as 
well  as  bearing  drainage  fluid  temperatures  were  also  measured  at  the  out¬ 
board  journal  and  thrust  bearings. 

Pump  inlet  pressures  were  measured  before  and  after  the  inlet  filter. 
Other  gages  were  installed  to  measure  pump  discharge  pressure,  interstage 
pressure  for  each  stage  and  bearing  lubrication  inlet  pressure. 

PUMP  TESTING  AND  RESULTS 

Four  Stage  Pump  Tests 

Testing  of  the  four  stage  hydraulic  pump  had  a  two-fold  purpose.  First 
to  prove  the  hydrodynamic  design,  and  secondly  to  demonstrate  an  endurance 
life  of  25  hours  with  a  pump  fluid  inlet  temperature  of  600°F. 

During  testing,  the  pump  achieved  a  peak  efficiency  of  62  percent.  This 
includes  all  mechanical  and  parasitic  losses  associated  with  the  pump  unit 
but  does  not  include  mechanical  or  thermodynamic  losses  of  the  air  turbine 
drive. 


After  the  final  thrust  bearing  configuration  was  established,  the  unit 
was  operated  for  a  total  twenty  hours  and  fifty  minutes.  Of  this  time,  fluid 
inlet  temperature  was  400°F  or  greater  for  fourteen  hours  and  twenty  minutes, 
and  600°F  for  six  hours  and  thirty  minutes. 


A  total  of  ten  short  test  runs  were  made,  and  were  terminated  for  a 
variety  of  reasons  associated  with  either  the  pump  test  or  the  pump  test  set 
up.  Table  I  lists  each  of  the  runs,  the  major  results,  and  reasons  for 
shutdown.  Test  data  for  each  of  the  runs  is  listed  in  Table  2  through  8  . 
Data  is  included  from  the  initial  break-in  runs  made  using  MIL-L-608IB,  Grade 
1005  lubricating  fluid  and  a  hydrodynamic  pocket  thrust  bearing. 

The  pump  was  completely  disassembled  following  Run  No.  3  and  Run  No.  7 
low  output  flows  were  observed  for  all  runs  subsequent  to  Run  No.  4,  and 
during  that  run,  flow  could  not  be  recorded  accurately  because  of  a  flow¬ 
meter  malfunction.  During  Run  No.  3,  when  full  output  flow  was  measured,  an 
overall  pump  efficiency  of  62  percent  was  achieved.  This  efficiency  was 
measured  as  follows: 

Turbine  shaft  power 


HP 


Wt6T 

TTT 


(  104. OH  195) 
177 


a 


I  15  horsepower 
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TABLE  I 


SUMMARY  OF  TEST  RUNS 


Test  Run 

Date 

Remarks 

1 

9-3-65 

Maximum  temperature  350°F,  flow  low 

2 

9-7-65 

Heat  Exchanger  Bypass  Valve  leak 

3 

9-0-65 

Rear  Seal  fai led 

4 

10-19-65 

Discharge  Valve  failed 

5 

10-21-65 

External  leak 

6 

10-25-65 

Low  output 

7 

10-27-65 

Shaft  seal  leak 

8 

12-1-65 

Seal  leak 

9 

12-2-65 

Test  stand  leak 

10 

12-3-65 

Fluid  breakdown 
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TEST  RUN  4,  5  AND  6,  OUTLINE  590020(L  208043) 
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TEST  RUN  1,  OUTLINE  580020  (L  2nf,043) 
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Estimated  from  previous  turbine  qualification  test 


TEST  RUN  10,  OUTLINE  580020(L  208043) 
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Pump  Hydraulic  Power 


HP  ,  Mjf  ,,  ,  65.9  horsepower 

♦Mechanical  Losses  =  6.0  horsepower 

This  efficiency  includes  mechanical  losses  at  the  pump  thrust  and 
journal  bearings,  the  nine  control led-gap  floating-carbon  bushings  and  the 
two  mechanical  face  seals.  It  also  includes  hydraulic  losses  attributable 
to  the  effect  of  the  bleeding  of  the  pump's  first  stage  to  provide  bearing 
lubricat ion. 

Since  no  attempt  was  made  to  determine  the  minimum  bearing  lubrication 
it  seems  reasonable  to  assume  that  efficiency  could  be  raised  even 
higher.  An  effort  was  made  to  get  some  feel  for  the  effect  of  bearing  flow. 
It  was  observed  that  a  slight  change  in  bearing  oil  flow  (first-stage  bleed) 
had  a  major  effect  on  first-stage  discharge  pressure  and  pump  total  dis¬ 
charge  pressure. 

Discussion  of  Test  Results 

In  discussions  of  tect  results  and  the  problems  encountered  during 
testing,  mechanical  design  of  the  equipment  must  be  considered  upper  most 
since  one  of  the  test  objectives  was  to  demonstrate  a  pump  life  of  25  hours 
at  high  temperature  operation  and  secondly  many  of  the  problems  and  delays 
were  directly  related  to  the  mechanical  design.  The  following  is  recapi¬ 
tulation  of  the  major  design  and  test  problems. 

I .  Thrust  Bearing 

The  initial  break-in  runs  were  halted  because  of  problems  with  the 
hydrodynamic  pocket  thrust  bearing  design.  Following  the  break-in  runs 
partial  disassembly  was  made,  and  it  was  found  that  both  sides  of  the  thrust 
runner  were  badly  worn.  The  mating  thrust  bearings  were  also  scored,  but 
th®  silver  plated  surface,  was  not  badly  worn.  The  inboard  journal  bearing 
plating  was  smeared  and  the  journal  showed  many  small  surface  cracks.  The 
journal  apparently  lost  radial  clearance  at  the  time  the  thrust  bearing 
overheated.  The  outboard  bearing  and  journal  was  in  excellent  condition. 

This  would  indicate  that  the  bearing  failure  was  not  caused  by  an  incompati¬ 
bility  of  materials.  Figure  78  shows  both  thrust  runner  surfaces,  the  thrust 
bearings,  and  both  journals. 


*From  past  data  recorded  during  many  development  and  qualification  tests  of 
AiFtesearch  Cabin  Turbocompressor  205900  which  is  the  cabin  turboccmpressor 
used  on  all  versions  of  the  Boeing  707  series  airplanes.  The  pump  drive 
turbine  is  the  turbine  section  of  this  well  tested  unit. 
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The  design  of  the  hydrodynamic  pocket  thrust  bearing  incorporated  a 
spherical  surtax  mating  with  a  spherical  seat  to  accommodate  slight  misa¬ 
lignments  and  surfaces  that  were  out  of  normal.  The  bearing  was  prevented 
from  rotating  by  a  single  loose-fit  dowel  pin.  From  the  almost  identical 
wear  patterns  on  both  sides  of  the  thrust  runner,  it  appeared  that  high¬ 
speed  operation  of  the  thrust  runner  caused  both  thrust  bearings  to  rotate 
about  the  retaining  dowel  pin.  This  resulted  in  loss  of  axial  clearance 
(nominally  0.006  in.)  as  the  bearings  moved  along  the  spherical  seat  and 
pinched  against  the  runner.  This  is  further  indicated  by  the  similar  wear 
patterns  on  the  thrust  bearings  at  identical  radial  positions  with  respect 
to  the  retaining  dowel  pin. 


This  particular  bearing  design  w a*  originally  selected  because  of  the 
low  losses  compared  to  ether  thrust  bearing  designs.  However,  rather  than 
expend  additional  effort  to  develop  this  rather  costly  bearing,  a  standard 
pivoted  shoe  type  thrust  bearing  was  installed  and  used  successfully  for  all 
subsequent  runs.  Figure  79  shows  this  bearing  following  Run  No.  3. 


2.  Carbon  Face  Shaft  Seals 


The  unit  disassembly  following  Run  No.  3  was  caused  by  failure  of  the 
outboard  shaft  seal.  The  seal  was  installed  to  duplicate  thrust  loads  that 
occurred  during  break-in  runs  when  the  outboard  bearing  lubrication  flow 
was  measured  at  ambient  pressure  rather  than  pump  inlet  pressure.  The  seal 
apparently  failed  because  of  high  temperatures  encountered  and  the  low 
lubricity  of  thr  MLO  8200  fluid. 


Local  heating,  as  a  result  of  the  seal  failure,  caused  burning  of  the 
outboard  journal  and  seal  runner.  Figure  60  shows  the  journal.  As  can  be 
seer  in  Figure  80,  there  are  additional  score  marks  in  the  area  of  the 
floating  carbon  seals.  This  was  apparently  the  result  of  journal  thermal 
expansion  caused  by  the  mechanical  seal  failure.  Wear  at  the  outermost 
bushing  was  measured  as  0,001  inch. 


During  the  test  runs,  the  need  of  a  outboard  mechanical  face  seal  was 
checked  by  venting  the  chamber  behind  the  seal  to  the  fluid  reservoir. 
Normally  this  chamber  was  vented  to  ambient  pressure;  no  change  in  thrust 
bearing  loading  was  noted  between  these  modes  of  operation.  Therefore,  it 
was  concluded  that  this  sea!  is  not  required  with  the  present  bearing 
design,  and  it  was  eliminated  for  ail  subsequent  runs. 

A  second  problem  with  carbon  face  shaft  seals  was  noted  following  Run 
No.  7  and  again  following  Run  No.  8.  The  answer  to  the  question  regarding 
the  effect  of  the  high  temperature  MLO  8200  fluid  interacting  with  ambient 
air  at  the  input  shaft  (see  Section  4,  Page  GS  was  noted.  During  Run  No.  7, 
the  unit  was  operated  with  a  fluid  inlet  temperature  of  600 °F  for  approxi¬ 
mately  two  hours.  A  shut  down  was  made  to  refill  the  CO 2  tanks  used  to 
provide  the  inert  atmosphere.  At  the  restart,  a  major  leak  was  noted  at 
the  input  sh^  ft.  Fluid  seepage  residue  had  built  up  to  the  level  of  the 
carbon  nose,  lifting  the  nose  from  the  seal  runner  and  permitting  free  flow 
of  the  fluid  past  the  seal.  Figure  81  shows  this  buildup. 
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81.  Carbon-Face  Shaft  Seal  After  Test  No.  7 


A  probe  was  installed  to  direct  inerting  gas  directly  into  this  area 
for  Run  No*  8,  however,  the  sludge  still  developed.  An  additional  high 
pressure  probe  was  then  installed  and  Runs  No.  9  and  No.  10  completed. 

3.  Low  Pump  Flow 

As  mentioned  previously,  for  most  of  the  test  runs,  pump  flow  was 
lower  than  design  flow  rate.  This  has  been  mainly  attributed  to  the  in¬ 
ability  of  properly  matching  impellers  and  diffusers.  Pump  tip  width  is 
only  0.040  inch.  The  unit  was  shimmed  during  assembly  to  match  each  of 
the  four  impeller  stages  with  the  mating  diffuser  and  to  provide  a  shroud 
clearance  of  0.019  to  0.021  inch  on  each  side  of  the  impeller.  However, 
because  of  the  method  of  shimming  (grinding  the  impeller  hub  to  suit)  and 
because  of  the  cumulative  effect  of  the  tolerance  stack  up,  it  was  ex¬ 
tremely  difficult  to  maintain  the  proper  clearances,  and  impossible  to 
check  for  tip  matching  and  shroud  clearances  after  the  final  assembly  ./as 
made.  In  this  respect,  it  is  suggested  that  future  designs  incorporate 
a  modular  principal  in  order  to  prevent  the  excessive  tolerance  stack  up. 

4.  Fluids 


The  final  run  was  halted  because  of  a  complete  breakdown  of  the  ML0 
8200  fluid.  In  addition  to  the  known  poor  lubricating  qualities  of  this 
fluid,  which  isamixtureof  alkyl  silicate  ester  compounds  has  several  well 
known  potential  problems.  These  are  thermal  degradation,  oxidation,  and 
susceptibility  to  hydrolysis.  The  effects  of  this  breakdown  are  an 
crease  in  acidity,  a  reduction  of  the  flashpoint  and  the  precipitation 
of  the  products  of  decomposition. 

The  mode  of  failure  either  during  or  following  the  final  test  run  wa. 
hydrol,sis,  noted  by  an  extremely  high  water  content  of  the  fluid  after  the 
test. 


As  a  result  of  the  fluid  breakdown,  both  journals  failed,  and  some 
damage  was  noted  at  the  thrust  bearing  runner.  Figures  82  and  83  show  the 
bearings  following  this  run. 

5.  Test  Setup 

Hany  delays  and  stoppages  were  caused  by  leaks  and  failures  of  test 
stand  components.  Although  these  problems  are  net  within  the  scope  of  this 
conflict,  the  area  is  mentioned  in  order  to  emphasize  the  need  ter  extreme 
care  in  design  and  selection  of  test  stand  components. 
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APPENDIX 


ORONITE  HIGH  TEMPERATURE  HYDRAULIC  FLUID  70 


Oronlte  High  Temperature  hydraulic  Fluid  70  is  a  new  silicate  ester  fluid 
developed  for  the  Oronite  Division,  California  Chemical  Company,  by  the 
California  Research  Corporation.  The  fluid  was  designed  as  a  replacement 
for  Oronite  8200  Fluid  in  the  B-70  Bomber  Program.  The  improvements  sought 
in  the  development  program  were  a  higher  operating  temperature  and  better 
shear  stability.  Because  of  the  stage  of  development  of  the  B-70  hydraulic 
system,  it  was  necessary  that  the  new  fluid  be  suitable  as  a  direct  re¬ 
placement  for  Oronite  8200  Fluid  in  the  system.  This  meant  the  new  fluid 
must  be  compatible  with  8200  Fluid  and  with  all  the  parts  in  the  system  and 
that  there  be  no  major  differences  in  the  physical  properties  of  the  two 
fluids.  All  of  these  goals  were  achieved  in  Oronite  70  Fluid,  and  addi¬ 
tional  improvements  were  found  in  other  70  Fluid  properties.  The  operating 
range  for  Oronite  70  Fluid  is  -65*F  to  +630“F  versus  a  range  of  -6 5*F  to 
+520°F  for  Oronite  8200  Fluid. 

The  attached  table  gives  comparative  test  data  for  Oronite  JO  Fluid,  Oronite 
3200  Fluid  and  MIL-H-5606  Fluid,  and  the  requirements  of  MIL-H-A446B  speci¬ 
fication.  In  the  following  discussion,  the  properties' of  Oronite  70  Fluid, 
Oronite  8200  Fluid  and  5606  Fluid  are  compared.  It  is  recognised  that  5606 
Fluid  is  not  a  high  temperature  hydraulic  fluid  and  it  is  not  designed  to 
operate  In  the  high  temperature  ranges  where  Oronite  70  Fluid  and  Oronite 
0200  Fluid  can  be  used.  The  comparison  has  been  Included  only  because  5 606 
Fluid  is  such  a  well  known  reference  in  the  industry. 

1.  Shear  Stability 

Oronite  JO  Fluid  possesses  excellent  shear  stability  as  measured  by 
the  sonic  oscillator  test.  Comparative  2-hour  sonic  oscillator  tests 
run  on  70  Fluid,  8200  Fluid  and  5606  Fluid,  70  Fluid  retained  98$  of 
its  original  viscosity  as  measured  at  210*F  versus  6j%  and  56$  reten¬ 
tion  for  8200  Fluid  and  5606  Fluid,  respectively. 

2.  Thermal  Stability 

Oronite  JO  Fluid  has  much  better  thermal  stability  than  8200  Fluid. 

70  Fluid  loses  only  13$  of  its  210*F  viscosity  when  held  at  600*F  for 
20  hours  versus  a  viscosity  loss  of  for  8200  Fluid  under  the  seise 
conditions.  Comparative  tests  were  not  run  on  5606  Fluid  because  its 
maximum  operating  temperature  is  only  275*F. 
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3.  Viscosity  Temperature  Properties 

As  shown  in  the  attached  table,  Oronite  70  Fluid  has  lower  Initial  over* 
all  viscosity  than  8200  Fluid.  However,  because  of  its  superior  thermal 
and  shear  stability,  the  viscosity  of  70  Fluid  holds  up  ouch  better  than 
the  viscosity  of  8200  Fluid  and  5^06  Fluid  after  a  relatively  short  per¬ 
iod  of  operation.  After  a  short  period  of  operation,  the  viscosity  of 
Oronite  70  Fluid  will  be  higher  than  that  of  8200  Fluid  or  5606  Fluid. 
Oronite  70  Fluid  does  not  contain  the  conventional  polymeric  viscosity 
index  improver,  and  for  this  reason,  the  fluid  has  better  thermal  and 
shear  stability  and  hence,  better  over-all  viscosity  properties. 

4.  Density 

The  density  of  70  Fluid  (0.953  gm  per  cc)  is  slightly  higher  than  that 
of  Oronite  8200  Fluid  (0.932  go  per  cc)  or  5606  Fluid  (0.860  go  per  cc). 
Weight,  of  course,  is  an  extremely  critical  factor  in  aircraft  design. 
However,  the  above  increase  of  2. 3$  in  the  density  of  70  Fluid  is  with¬ 
in  acceptable  limits  for  the  B-7C  Program. 

5.  Flash  Point  and  Autogenous  Ignition  Temperature 

The  flash  point  of  70  Fluid  is  43Q*F  versus  395°F  for  8200  Fluid  and 
225“F  for  5606  Fluid  The  autogenous  ignition  temperature  of  70  Fluid 
is  735#F  versus  760°F  for  8200  Fluid  and  475*F  for  5606  Fluid. 

6.  Hydrolytic  Stability 

Oronite  70  Fluid  and  8200  Fluid  both  meet  the  requirements  of  MIL-H- 
8446b  hydrolytic  stability  test?.  The  appearance  and  weight  charge  of 
the  copper  is  essentially  the  fame  for  both  fluids.  However,  the  vis¬ 
cosity  increase  and  acid  number  of  the  fluid  is  slightly  higher  for  70 
Fluid  than  for  8200  Fluid.  No  comparison  is  made  here  with  5606  Fluid 
since  the  hydrolytic  stability  tests  of  a  petroleum  oil  is  meaningless. 

7«  Oxidation  and  Corrosion  Stability 

Oronite  70  Fluid  glues  better  results  than  8200  Fluid  on  the  MIL-H-6446B 
oxidation  and  corrosion  stability  test.  The  weight  change  and  appear¬ 
ance  of  the  metal  specimens  is  essentially  the  same  for  both  fluids. 
However,  the  viscosity  change  and  final  acid  'number  of  the  fluids  is 
better  for  70  Fluid  than  for  8200  Fluid.  The  viscosity  change  of  70 
Fluid  is  only  +10$  versus  -19$  for  8200  Fluid  and  the  final  acid  nusber 
of  70  Fluid  is  0.1  versus  0.6  for  8200  Fluid.  Direct  comparisons  with 
the  oxidation  corrosion  stability  of  5606  Fluid  cannot  be  made  since 
this  fluid  is  normally  tested  under  less  severe  conditions.  In  the  MXL- 
H-3446B  oxidation  and  corrosion  test  used  on  70  Fluid  and  8200  Fluid, 
the  sample  is  held  at  400*F  for  72  hours.  5&06  Fluid  Is  normally  tested 
at  250*F  for  168  hours.  Under  these  latter  conditions,  5606  Fluid  shows 
a  viscosity  change  of  +3*4$  and  a  final  acid  cumber  of  0.1$. 
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8.  Vapor  Pressure 

. .  ♦ 

The  vapor  pressure  of  70  Fluid  is  markedly  better  than  8200  Fluid  or 
5606  Fluid.  At  k00*F,  the  vapor  pressure  In  millimeters  of  mercury 
for  70  Fluid  is  0.8,  for  8200  Fluid  -  3-0,  and  for  5606  Fluid  -  220. 


9.  Lubricity 

Oronite  70  Fluid  gives  essentially  the  same  results  on  the  four  ball 
wear  test  as  Oronite  8200  Fluid  indicating  similar  lubricity  properties 
of  these  tvo  fluids.  Four  ball  wear  test  conducted  at  167*F  showed 
that  5606  has  slightly  better  lubricity  than  70  Fluid.  However,  four 
ball  wear  test  cannot  be  conducted  on  5606  Fluid  at  h00*F  since  this 
fluid  will  decompose  at  that  temperature. 


10.  Foaming  Tendency 

One  of  the  outstanding  characteristics  of  Oronite  70  Fluid  is  that  it 
does  not  form  any  foam.  In  all  of  the  foam  tests  which  have  been  run 
on  this  fluid,  the  volume  of  foam  formed  is  too'samll  to  be  measured. 

The  foam  breaks  immediately.  In  comparative  foam  tests,  Oronite  JO 
Fluid  gave  no  foam,  8200  Fluid  gave  350  millimeters  of  foam  and  5606 
gave  U5  millimeters  of  foam. 

In  summary,  Oronite  High  Temperature  Hydraulic  Fluid  70  has  several  advantages 
over  Oronite  High  Temperature  Hydraulic  Fluid  8200.  It  has  a  wider  operating 
temperature,  better  thermal  acidative  and  6bear  stability  than  8200. Fluid, 
and  is  a  nonfoaming  fluid.  Oronite  70  Fluid  has  better  over-all  viscosity 
properties  than  8200  Fluid  and  has  a  lower  vapor  pressure  and  higher  flash 
point  than  8200  Fluid.  The  minor  disadvantages  jf  slightly  poorer  hydrolytic 
stability  and  higher  density  are  greatly  offset  by  this  new  fluid's  many 
advantages . 

Oronite  70  Fluid  is  a  new  concept  in  high  temperature  hydraulic  fluid  design. 
The  absence  of  a  conventional  polymeric  viscosity  index  improver  greatly 
enchances  the  fluid's  properties.  The  conventional  polymeric  viscosity 
index  improver  is  the  weakest  link  in  a  hydraulic  fluid  from  a  thermal  and 
shear  standpoint.  Because  of  the  absence  of  the  viscosity  index  improver 
in  70  Fluid,  this  fluid's  viscosity  properties  are  not  degraded  over  a 
period  of  time  by  thermal  or  shear  action. 
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PROPERTIES  OF  TYPICAL  ORONITE  HIGH  TEMPERATURE  HYDRAULIC  FLUIDS  70  AND  8200 
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COMPARISON  OF  PROPERTY  REQUIREMENTS 
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THERMAL  STABILITY  OOMPARISOH 
OROr.TTE  HIGH  TEMPERATURE  HYDRAULIC  FLUIDS  70  AID  8200 


Percent  Viscosity  Change  —  Measured  at  210*? 


70  Fluid _ 8200  Fluid 


Temperature 

and 

Tine 

Fluid  as 
Removed  from 
Chpsuie 

Fluid  After 

Removing 

Volatiles 

Fluid  as 
Removed 
frc-i  Capsule 

"  iKdrm£ — 

Removing 

Volatiles 

6C0*F 

6  hours 

1 

-7.1 

♦6. 5 

-50 

20  hours 

-li 

+2.6 

.64 

-59 

40  Lours 

-32 

-10 

-72 

62 

630*F 

*4.6 

6  hours 

-13 

-63 

-54 

20  hours 

-38 

-5.6 

-73 

mr.7 

s  ; 

40  hours 

_ 

-49 

+26 

mm  •• 

m  m>  m 

In  the  above  test  10  ml  camples  of  the  hydraulic  fluiu.  were  sealed  In  stainless 
steel  capsules  under  an  atmosphere  of  helium.  The  fluids  were  heated  to  600*? 
and  630 *F  for  the  tine  periods  indicated  above.  TLe  samples  were  removed  from 
the  capsuJ  ?s  and  the  viscosities  measured  at  210*F.  T  e  volatiles  formed  during 
the  high  temperature  period  in  the  capsules  were  remo*  d  by  evaporation  and  the 
viscos  '  ties  -ere  remeasured.  The  results  reported  ibove  are  the  per  cent 
changes  in  the  viscosities  from  the  original  viscos ttlee  as  measured  at  210*F. 
The  original  2iO*F  viscosities  were  for  70  Fluid  7.56  centistohes  and  for 
8200  Fluid  11.27  centistcJcea. 
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FIGURE  84 


PROPERTIES  OF  ORONITE 
HIGH  TEMPERATURE  HYDRAULIC 
(and  reference  fluids 
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FIGURE  84  (continued) 


COMPUTER  PROGRAMS  FOR  PERFORMANCE  CALCULATION  (PROGRAMS  A-0545  AND  A-CS46J 

The  following  paragraphs  present  the  computer  program  w^stter,  for  the 
method  of  performance  celeulft ion  developed  previously. 

Program  A-0545  calculates  pumc  designs  for  all  combinations  of  inps:t 
values  for  N,  ,  p.,  and  If  desired,  it  will  then  ca Icy! ate  off-design 
points  for  different  values  of  MQ  and  M/N^  Automatic  plotting  of  AP  and  Tj 

curves  for  pump  off-design  is  also  possible.  Programs  A-C546  and  D-0100  are 
chained  to  A-0545  for  this  purpose. 

Input  for  A-0545 

If  only  the  pump  design  point  program  is  used,  !l  cards  of  input  are 
needed  for  each  case.  If  the  pump  off-design  program  is  used.  14  cards  of 
input  are  needed  for  aach  cate.  Input  for  pump  design  point  program; 


Card 

No.  1 

format 

72H 

A 1 phanumer i c 
"i"  in  column 

heading  card 

1 

Card 

No.  2 

FORMAT 

5FI0.5 

Y0»  !  b/cu  ■'{. 

Dens  i  t  y 

V  ,  C‘j  ft /sec 
o 

,  F  I  ow 

APo.  psi, 

Pressure  rise 

Vtt,  ft2/ sec , 

K i nemaf i c  v i scos i t y 

Test 

-  0  if  pump  design  point  program 

i  if  pump  off -design  point  program 

Card 

No .  3 ,  Inlet 

FORMAT  110 

.  tfio.5 

Case  No 

1  NPSH  Given  -  none av <  tai i ng 

2  0.  given 

i 

3  Best  performance 

4  Diame  ers  and  NPSH  given  -  cavitating 

operat i on 

°H ’  in‘ 

Hub  diameter 

NPSH,  ft 

Net  positive  suction  head 

0j .  deg 

Inlet  flow  ang 1 e 
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K. 

i 

D . ,  in. 

i 

A,  I/ft 

f 

K 

4 

4.  Card  No.  4 

c . 

I 

5.  Card  No.  5,  N  values 

No.  of  N  values 
N,  RPM 
N? 

N3 

etc. 

6.  Card  No.  6,  Impeller 

?' 

?* 

?3 

?4 

7.  Card  No.  7,  values 

No.  of  g?  values 

0i(») 

MO 

3 ) 

etc . 


Hub-blockage  K.  - 

Inlet  tip  diameter 

Density  function  constants 

Blockage  factor 

(max  10  values) 


Loss  constants 


Impeller  tip  vane  angle 


FORMAT  FI 0.5 


FORMAT  110,  7FI0.5/ 
8FI0.5 


FORMAT  4FI0.5 


FORMAT  110,  7FI0.5/ 
8FI0.5 
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8 .  Card  No .  8 ,  p  va  I  ties 

No.  of  p  values 

M-'  . 

U? 

ULJ 
etc . 

9.  Cord  No.  9,  \  values 

No.  of  X  values 

X, 

X* 

Xj 

etc. 

(0.  Card  No.  10,  Diffuser 
Case  No. 


Y  >  deg 
?* 


?7 

Dj 


Slip  factor 


Flow  factor 


JUj® 


juj 

u2 


FORMAT  HO,  7FI0.5/ 
8FI0.5 


FORMAT  HO,  7FI0.5/ 
8FI0.5 


FORMAT  HO,  6FI0.5 


1  -  Vaned  di f fuser 

2  -  Scrol I  di f fuser 

3  =  Coll ector  r i ng 

4  -  Vaneless  diffuser 
Throat  s i de-to-wi dth  ratio 
Wall  angle 

Boundary  layer  factor 
Boundary  layer  factor 
Loss  coefficient 

Vaneless  space  diameter  ratio 
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II.  Card  No,  II,  Outside  Losses 


FORMAT  2FI0.5 


Wheel  friction  loss  factor 
Leakage  loss  factor 

Input  for  Pump  off-design  program: 

Cards  I  thru  II  of  previous  section  are  needed  in  addition  to  the  following 
12.  Card  No.  12  FORMAT  8FI0.5 


a2,  I/Deg 
as,  I/Deg * 

*c 

?! 

?3 

13.  Card  No.  13 

u 

?• 

v*  -ft1/ sec 
y  -lb/cu  ft 
NPSH,  ft 


Incidence  loss  coefficients 

Cavitation  loss  coefficient 


FORMAT  5FI0.5 

Diffuser  incidence  loss  coefficient 

Kinematic  viscosity  for  part  load 
Density  for  part  load 

Net  positive  suction  head  for  part  load 
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14.  Card  No.  14 


FORMAT  6F10.5 


M  Mass  flow  ratio 

o 

M 

o 

(Mo^N 

Speed  ratio 

A(N/Nq) 

Input  for  A- 0546 

If  plotting  of  pump  off-design  points  is  desired,  two  blank  cards  plus  the 
following  cards  should  be  inserted  behind  the  data  deck. 

1.  Card  No.  I  FORMAT  3110 

No.  of  problems  to  plot 

No.  of  N/Nq  values  for  1st  problem 

No.  of  Mq  values  for  1st  problem 

2.  Card  No.  2  FORMAT  2110 

No.  of  N/Nq  values  for  2nd  problem 
No.  of  M^  values  for  2nd  problem 

3.  Repeat  Card  No.  2  for  as  many  problems  as  are  to  be  plotted. 

Output 

The  output  consists  of  three  parts  (a  copy  of  a  sample  output  is  attached)^ 

1.  The  input  is  printed  out 

2.  For  a  particular  N,  P*,  \k  and  X  the  following  quantities  are  printed  out: 
a.  Inlet  Case 


N,  RPM 

O. ,  inches 

Pi»  deg 
NPSH,  ft 


Speed 

Inlet  diameter 

Inlet  flow  angle 

Net  positive  suction  head 
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V,,  F PS 

Inlet  relative  velocity  inducer 

Suet  ion  dens i ty  rat lo 

Impeller  Input 

fii>  Deg 

Impeller  tip  vane  angle 

51  ip  factor  (  Cu?/CiJioo) 

Flow  factor  (  Cm,/u») 

Impeller  Answers 

1  q 

Head  factor 

Of,  inch 

Impeller  tip  diameter 

U i,  FPS 

Impeller  tip  velocity 

?. 

Impeller  relative  velocity  ratio 

z, 

Number  of  vanes 

2  B, 

Impeller  width  ratio  (=  b*/D») 

R  . 
ei 

Impeller  channel  ave.  Reynolds  No. 

AqSF 

Impeller  friction  loss 

^OL 

Diffusion  loss 

AqRC 

Recirculation  loss 

T'imp 

Impeller  efficiency  (Internal) 

Diffuser  Answers 

If  Case  No.  1  -  the 

following  line  is  printed: 

1 •  Case  No. 

Oiff"Ser  inlet  Reynolds  No. 

«e2 

Impeller  tip  Reynolds  No- 

Zi 

fl>>».,'er  of  diffuser  vanes 

*Sc 

Scrol 1  frict ion 

*dif.  2 

Straight  diffuser  loss 

Total  d  ffuser  loss 
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Application  of  the  IBM  Programs  to  Performance  Calculation 


1 .  Generalized  Performance  of  Centrifugal  Pumps 

In  order  to  compare  the  calculation  results  with  statistical  average 
and  predicted  optimum  performances,  a  number  of  data  were  gathered  on  high- 
Reynolds-number  pump  performances  and  plotted  in  Figure  42.  For  comparative 
purpose,  three  points  were  calculated  by  the  present  computer  program  using 
water  as  fluid,  assuming  no  cavitation  and  a  Reynolds  number  close  to  I07. 

The  specific  speeds  of  the  calculated  samples  were  500  to  1000  and  2000.  The 
results  of  the  calculations  were  also  plotted  in  the  diagram.  The  points  fit 
the  top  efficiency  curve  quite  well;  the  program  therefore  might  be  used  to 
predict  the  performance  of  advanced  pump  designs  with  confidence. 

2 .  Performance  Characteristics  of  the  LH?  Pump 

The  performance  calculation  procedure  developed  previously  has  been 
applied  to  establish  the  characteristic  curves  of  the  pump  designed  and  built 
in  this  program.  Figures  38  and  39  present  the  characteristics  over  a  wide 
range  of  speed  and  flow. 

Execution  Instructions 


Program  A-0545  consists  of  three  chains.  Chains  2  and  3  are  Programs 
A-0546  and  0-0100  and  are  used  only  for  plotting. 

Tape  22  is  always  required  by  Program  A-0545.  If  plotting  is  called  for, 
tape  22  will  be  the  plotter  tape.  When  plotting  is  done,  tapes  21  and  II  are 
used  as  work  tapes. 

Each  of  the  programs  just  mentioned  takes  approximately  1-1/2  min  per 
case  with  150  solutions. 

Following  this  page  are  given: 

1.  Flow  chart  of  program 

2.  Fortran  listing  of  program 

3.  Sample  input  data  listing 

4.  Sample  output  listing 

5.  Samp  I e  plots 
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Figure  85.  Flow  Chart  of  Program 


Jim  AG566  pump  DESIGN  plot  program  with  garlic  chain 

123  FORMAT ( 5F 14 . 7 ) 

126  FORMAT  I  ItiMP’iO.O  10.0  ,  I  2  f  IM*  ,  4X,  1H1 , 4X,F4. 2  ,  6X  ,4F1 0.  1 ) 

127  FORMAMloX,  2F  10.71 

128  FORMAT ( 2 HH  ,2F7.1) 

129  FORMAT  I  1  UIX  0.0  .0,20X,18H  MOIN/NOl) 

130  FURMATI  3CHAC-5.  9.0  .12  PUMP  -  PART  LOAD  PERFORMANCE/! 

131  FURMATU3HY  0.0  0 . 0 , 2 J X , 1 8H0VFR ALL  EFFICIENCY) 

332  FORMAT!  13HY  0.0  0 . 3 , 20X ,  18HPKE  S  SURE  RISE,  PSD 

l  33  FURMATI  1  2ltG  2.0  9.31 

136  FORMAT (  WHAOl .5  R.6  .12  FLOW,  CUFT/SFC  «  F7.3/I 

136  FURMA  T I  IKE 1 

137  FOHMAT I 10X, F10. 7.F10.4) 

1  38  FURMATI  31  10 1 

DIMENSION  ANSI  I  2  I , YOR I  2  1 , OE L T Y ( 2  1 

X  UR  -0.0 

DEI Tx=2.0 

YOR 111=0.0 

YURI21O.0 

OEL  TYI  l  )  =  l .0 

OtL  TV! 2 1=6000. 

X UR  IC»N  =  O.C 


YUR IGN=0.0 
REWIND  65 

R  t  AO  INPUT  TAPE  6  1 ,  13R,K1,K2,K3 
WRITE  nuTPUT  TAPE  45,I28,X0RIGN,Y0RIGN 
00  HI  6  1=>,KI 
IFIL-U  779,  779,  778 

778  READ  INPUT  TAPF  41 ,  1  38,  K2,K  3 

779  00  313  K=l,2 

GO  TO  I  784,  7P0, 7HC, 780,  780,  780)  ,L 

780  GO  TO  I  785, 7R1 1 ,K 

781  K  4  =  K  2  »K  3 

DU  7B2  I  =  l  *  K4 
7R2  BACKSPACE  69 
GU  TIT  785 
786  REWIND  69 
785  00  HO 6  I = 1 , K2 
RCWINO  67 
J  =  0 

DO  3C3  I  *  !. ,  K  3 

READ  INPUT  T AIT  68, l 25, VO, ANSI! 2 ), ANSI  I l ) ,ANS, ANNO 
IF { ANSI  IK  1 } BC  3,  802,  802 

802  J=J»1 

WRITE  OUTPUT  TAPE  *.  7 , 1  2  5  ,  ANS  ,  ANS  1 1  X  1 

803  CONTINUE 

WHITE  OUTPUT  TAPE  6  5 , 1  2  6 ,  J ,  ANNO  ,  XTTR  ,  OE  L  T  X ,  YOR  I  K  1  ,  OE  L  T  Y  I  K  ) 

0  EW  I, NO  67 
00  H06  1  =  1 ,  J 

R  E  AO  INPUT  TAPF  6 7 , l 25 , ANS , ANS 1 1 K ) 

GO  TO  (806,805)  , K 

806  WRITE  OUTPUT  *APE  6 5 , 1 27 , ANS , ANSI! K 1 
GO  TU  806 

805  WRITE  OUTPUT  TAPE  6  3 , 1  3  7  ,  AN  S  ,  AN  S 1 1  K  )  v- 

806  CONTINUE 

WRITE  OUTPUT  TAPE  6  5,129 

WHITE  OUTPUT  TAPE  65,1  30  ' 

jU  Tf)  1808, 809  1  ,K  vy-*"' 


\ 
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808  WHITE  OUTPUT  TAPf  45, 1U 
GO  TO  810 

809  WHITE  OUTPUT  TAPE 

810  WHITE  OUTPUT  lAPt  45,131, 

WHITE  OUTPUT  TAPE  45,134,V9 

XQK  IG'1*  18  . 

B  15  WHITE  .HI  T  POT  TAPF.  4  5 , 1  ?  8  ,  XO»  !GN ,  VOR  I  ON 
816  CONTINUE 

WHITE  OUTPUT  TAPf.  45,136 
ENOFILE  4^ 

REWIND  45 
REWIND  46 
CALL  CHAIN! 03,461 
END 


AO  5  Ah 
AO  54  6 

AD  5  A  6 


•  JljM  A 0  V<  !>  Kl  I.PT  /MX  0  325/, 4  1*1  JM I*  PAR  I  LOAD  PERFORMANCE 

•  t  CUMP  A0V*‘> 

C  G  t  N  H<  A l  PUMP  f  F  F  |  C  I  F'Jt.Y 

DIM'  ‘I,  IHM  AlAMDAI  10  > ,  A4UII1  ■  ),fUTA7I(I0)  ,ANIT  1C) 

100  FOKMAI  {  I  10 1  If  1  0. 5/  C  RF  1 0.‘» )  ) 

1DI  FORMAT  I  nr  10. S) 

10?  FORMAl  (  13, 6f  10.  3,  10X,  3SIIINLE T  CASE  NO.  .  N  ,  0!  ,  IH  ,  NPSM ,  W  J  ,  L  E  F  A  ) 

103  FORMAT  (5  TOO  A  SCROLL  Of KL IKON,  04/03,  OGE  S  NOT  CONVERGE,  D4/DT-1 
1  *FH. 1,110,  D4/DJ-2  *r«.l,l?M,  ALPHA-3  =  EB.5) 

104  FORMAT  ( 36H0  CASC  4  Ci  DOFS  NOT  CONVERGE  Cll*E10.5,6H,  CIP’TIO.S 

1) 

105  FORMAT  (2?H0  INPUT-  DFNS  I  TV  *F7.3,f)H,  FLOW  *F7.3,15M,  PRESS.  R 
1ISC  =16.1,1511,  KJNtM.VISC.  =F9.6) 

107  FORMAT |4:JN0  IMPTLLFR  /I  IS  GRFATFR  THAN  100,  /I  =riO.S,6U,  Q  *  FI 
10.5,7M,  D?  *  F10.5.9H,  /FTAl  =F10.5) 

108  FORMAT  (Fh.1,?F 1Q.3,4uX, 77HIMPFLLFR-  KETA-2,  MU,  L AMOA/ F h . 2 , F 1 0 . 2 , 
IF  10. 1,  F1C.4  ,F  ID,  l  , 2  OX , ?uMU,  02,  U?  ,  ZFTA1,  Z  1  /F  3 . 5 ,  F  10.  1 , 4F  l  0 . 5 , 1  0 
2X 34002,  RFt,  OOSI  ,  0001. ,  OORC,  CTA-!mP) 

109  FORMAT  ( 4  7H0  DIFFUSLR  L?  NOT  OETwEEN  6  AND  16,  ZETA  01  *  F8.),12H 
1,  ZFTA  02  *  Fli.l,7H,  12  =  Ffi.2) 

110  FORMAT ( ?F 10.5, I 3H  0403,  D4D3N) 

114  FORMAT  (  3  4  HO  HTFIT  OOFS  NOT  CDNVFRGF  ,  HTOT  *  F10.5,  11M,  MTOT?  =■  F 
110.5) 

J.  1 5  FORMAT  (  18  ,?F  10.0,  T  10.2,  IF  10. 5,  7X45HD1  miSER  CASE  NO.  ,  R  F82  ,  Rt  2  ,  Z  2  , 
l  dus  c,  non  if  ?,  non » 

116  FORMAT  118, 4F 10. S, 30 X , 4QHD I FF USER  CASE  NO.,  D023,  n034,  D045,  000) 

117  FURMAT  I8X,F10.5,F10.3,  3F 10. 5, 2CX, 2 9MF T A-MYO ,  HP",  DUH  E ,  DQL,  ETA 
1) 

118  FORMAT  (72H 

1  ) 

119  FORMAT  (?9il0  0  DUES  NOT  CONVERGE,  01  *  F10.5,7H,  02  =  F10.5.7U,  MU 

1  =>  r;o.  5) 

120  F0RMATI5F 14.7) 

121  F0RMAT(F8.2,2F10.4,FIC.5,FI0.?,10X26H  MO,  0,  DFLTA-R,  ETAH,  HO/fiX 
l,F10.3,FlC.5,ri 0. I, 30X, 16M0ELTA-P,  ETA,  HP) 

122  FORMAT  I6X, 7MN/N0  *  F8.3) 

123  FORMAT ( 24H1  PART  LOAD  PERFORMANCE/) 

124  FORMAT  (36H0PART  LOAD  CIP  DOES  NOT  CONVFRG  C  1 1  =  F  1  0. 5 ,6H ,  CI2=ri0.5 

1) 

C  INLET 

REWIND  48 

1  READ  INPUT  TAPE  41,118 
WKI TE  OUTPUT  TAPE  42,118 

C  TEST  =  0  IF  PUMP  DESIGN  POINT  PERT,  =  1  .  IF  PUMP  OTF  DESIGN  PCRF 

READ  INPUT  TAPE  41,101,  GAMD, VO,DPO , VS  F  AM , TEST 

I F( GAMO«VO*DPO) 201 , 2GC, 201 

200  ENDFUE  48 
REWIND  4 H 

CALL  CM A  I N I  02 ,46) 

201  REAn  INPUT  TAPE  4 1 , 1  DO , I C A  ST , DM,  ANP SH , HI , AK  i  , DI , A , AK , C I 

READ  INPUT  TAPE  4  1 , 1 00 , N1 , I  AN  I ( I ) , I  *  1 , N1 ) 

READ  INPUT  TAPE  41,101,  X 1 , X? , X3 , X4 

READ  INPUT  TAPE  4 1 , l 00 , N? , ( nr T A 2 [( I  )  ,  1  =  1  , N? ) 

READ  INPUT  TAPF  4  1 ,  1 00 , N3 , ( AMU  II  I ) ,  I  =  1 , N3 ) 

READ  |N»UI  IAPF  '•!,  1  00,  N,  J  A  L  AMD  A  (  I  )  ,  I  =  1  ,  N  ) 

Rt AD  INPUT  TAPE  4 1 , 1 OC , I C AST  2 , 7 F T AD , G AMMA , X5 , X6 , X7 , D3D? 

SfcAD  INPUT  TAPT  4l,10l,X8,X9 
IF  I  TEST)  3, 3,2 
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2  READ  INPUT  TAPP  41,  t01,XlPl,M,A2,A3,PSIC.X2PL,X3PL,X4Pi.,PS!5,X8)*l. 

l,X9Pl  ,V$T  APL .DAMP, ANPSP  *054' 

f»£4D  INPUT  TAPE  41  ,  101 .  AMQO.'JELTM,  AMN,  ANNOO  ,DFL  TAN,  4NN0N 

3  G  AMK  =  F,  AXM  A*  .0  17453? 0 

TG AM- S I NF  (  GAWK  )  /COSE  IGA, MM  1 
CUNS4=2./  TGA1U.5 

4  rfKITE  OUTPUT  T  API  <,  2  ,  1 0  $ .  GAMO ,  VO  ,  OPO  ,  US  T  Aft 
DU  10  K 1  -  1 , N 1 

AN=AN1 ( K 1 ) 

GU  TU  (5, 6,0, 9)  .ICASE 

5  CONSl=(  12.1  75  •  A  \|P$H  )  ••  .  5 

01*1  *>?<».  *  VO/ (  T .  141 59?7»CF)\IS  l  1  *!)il»OH  )  ••  •  5 

8  I  -  AT  ANF  (3.14  11  )2  7»IJ|»AN/<  7  20.*C()NS1  i  > 

rfi»(  3?.  1  75»ANP$ti.  (  AN.  J.  141 1 ')?  7«  1)1  /7?0.  i  ••?)•<■.  1 

8  2  *  H  I  •  S  7 .  T  ^  S  7  ft 

GO  TU  lr> 

6  8  J  R  *1)  I  »  .0  1  745  3?‘> 

Ci»I=r.osr(  niRi 
CTHI=CHI/SINI  (A  IP/ 

7  01  =  144354. 36*  V  P/I  CTllI*  AN  »AKI  ))•».  33333 
Jl=3. 1415927* AN *01/720. 

C  I  *U1  •CTI*.  1 

wl*ui/siori  itir) 

ANPSH  =C I >01/64.34 
GU  TO  11 

e  81*60. 

GU  TU  6 

9  AKJ*l.-|OH/f)|)* Oh/D  1 
JI»0!*AN» 3. 1415927/720. 

CONS  1  =  144. •VO/IFI«(,M«OI  »AK  I  •  .70  5  390  2) 

J-l 

CU  *  CONSl 

IF  I CCNSl*C0NSl/64. 3 5- ANPSH) 14,  14,10 

10  C I  2  =  CONS  1 • ( l . ♦ A  » I C I  1 • C I  1/64, 3 5- ANPSH) )  •  •  AX 
|F<  A0SMCI1/C  1 2-1.) -.000 1)1  4,  14, 11 

11  J*J*1 

IF  (.1-50)  ’2.12,1  3 

i?  cn=ci2 

GO  TO  10 

13  WRITS  OUTPUT  TAPE  42  « 104 <  C  1  1  ,C  I  2 
CI1-ICII*C12)/?. 

14  81  =  6  ANFUU/CI1  i 
W  I  *li  I  /  S  I  NF  I  B I  ) 

ZETA*(1.+A»(CI1 »Cf 1/64. 35- ANPSH) }»«AK 
0 i  =  B I *5 7 .295  78 

15  WHITE  OUTPUT  TAPE  42,102, ICASE, AN, 0! ,81 , ANPSH, Ml ,  Z£  T  A 
1MPEUEK  CALCULATIONS 

B  l  R  =  l)  I  •.0174637  9 
C5l«Ct)SF(Hl#i 
T(31=SINF|!3!H)/tl3I 
00  50  K  7 *  I , N 2 
8ETA2=tlFTA2l  1K2) 

B  T 2  =Bh T  A2 " . 0 1 74 5  3  29 
$82  =  S  I NF  I  itT  2 ) 

CB2=CU$F ( BT2 I 
f  02  - SB2/C02 
OU  50  X  3  = 1 ; N  1 


AMU* AMU  I I  K  1  ) 

CUNS16-XU  -1  .5  70  r  )63-ArtU»CH2/I  l  .  -  AMU  ) 

DU  50  I  =  l ,  N 

CUNSU*  ALATH'Ai  j  )  »Al  AMDAl  I  ) 

U-AMU* ( 1 .-ALA^HAJ  I  )»I02  ) 

CUNS  ;» CCOUS*' ♦0»0)/(2.«0* 

CUNS‘»*  I  CUMS»i ♦0»U)  •  •  . 

C JN$6“  AL  AMI)A{  I  )  /  CONS'! 

C0NS1  /* I CU:4St' ♦  I  l.-O )••?)**. 5 
CQNS18*! * -CON SI 
J2-  1 

Hi)*  1  44  .  »l)PO/GAMi» 

E  T  AHYD* . 9 
HTOT-HO/E  r AHY 0 
16  U2*  I  32.1  75-MH1T /())  --.5 
02=  720. *U2/(  1.14150  2  7  •  A  N  ) 

W2*U?»C0NS1 / 

V2*D!/n? 

ZETAl=W?/vU 

Zl«CUN$16»t l.*.Otf/ZETAl )/( 1.-V2) 

IFIZl-100.  )2l,?l, ?0 

20  WHITE  OUTPUT  T  A  HC  42  » 10  7, Z 1 ,U,D2  .ZETAl 
GU  TO  50 

21  DEL  1  * . 8» I  02 -OH) M  CO  1 ♦CH 2 ) 

82  -144.  *  VO  /  (  .  0  •  AL  AMDAl  I  ) -UP- (0202  •  3.  1  41  50-02*  „  03*  Z  1  /C  82  )  ) 

0HY01  =  T)2  •  < 1  ./(/  1/(3.  1415927M»?)*l./Cl)2MV2/t  2./U.-V?)*  1 2. •/ 1/1 3.1 
141592 7 »ll  .*V?  i  ))•{  l  .*Tft  I»THI  -1 1  .*V2»V2)/2.  )**.5>  I 
WAVE*I(Wl»Wl»W2-W?)». !>)••. 5 
REl=OHVm-WAVE/l  l  2. -VST  AH) 

F  1  =  .U462/HC  I  »•.  2 

OQSF=X?»OEL 1 -WAVC^WAVE-Fl/I DHYD 1»U2»U2«0) 

DELTA=l.-/ETAl  +  U»(Cni+CU2)/I2.»WI»I  Zi « 11 .- V2 ) / 3 . 1 41 5 927  ♦2.*V2)/U2 

1) 

DUnL  =  X?»l DEL  TA»WI/U2)»»2/Q 
DQRC*X4»DELTA»0/ALAMDA( I) 

El  A  IMP -ICON  SI .  A-DOOL ) /ICONSlA+OQSFtDORC ) 

01FcUSEl< 

GO  TO  *  22,  37  ),  I  CASE  2 

VANE!)  DIFFUSER  W/OONICAl  SIDEWALLS 

22  RE2a02»U2/l  12.-VSTAK) 

REfs2=2.*B2*HF:  2*  CONS  5 
F  2* .0462 /RE  H? .  2 
DcL  * .  1 

J  l  *0 
J*  1 

ZETAU1-/5TA0 

23  CONS  7  -  /  F  TADl*TC,AM 

RR.O *  1  .04*  . 693  33 *  / F  T  A01  •  02  »  (  3 «  ♦?  .  »f.ONS  7 )  / 1  2 .  ♦CONS  7 ) 

FH-CONSA  ♦F?»IZETAD1*I2.*  T  G  A  N )  *• .  5  •  L  OG  F  M  l  6  T  40 1 -CB2  ♦  .0  2  ) 

1.02)  "C0NS4-LOGF  (IH2»IZFTADl»TGAM»2.  )  * .  02*  TGAM)  /  ( 2 « •  8  2  +  .  G  2  *  T  G  A  M  ) ) ) 
CM2Cli*t!R0  »F II*  CONS 5/0 

Z  2*  6. 2!)  3 1  85  3 -CM 200/  I  / 1  T  ADI  *  H2*  l  2.  ♦Cl,Ni7)  ) 

J  1  -  i  i  ♦  1 

IF  1  J1-50J  24,24,  34 

24  IF i / 2 - 6 . 0 ) 2 6 ,  <6,25 

25  IFI Z2-1A. 01 3< , 36,29 

26  GO  TO  (2H,?H,2?),J 


27  DLL  -  HI  I.  •() .  I 
?fl  J  = 

Z  l.  T  A i .  |  =  / 1  I  Al>  1  -  1)1  L 

Go  to  ?3 

29  GO  10  (  31  ,  3  41  )  ,  J 

30  OIL  -  Of  L  *1).  I 

31  J  *  4 

n. t  Ant  =  /m  aim  ♦  dm. 

GO  fO  2  4 

34  W  K  I  It:  OH  I  PUT  f  AIM  4  ?  ,  I  O') ,  l  M  ,0),  /  t  T  AO  1  ,  l  ? 

Go  ro  so 

U  DfcL?M>2»(  .'  //CONS/.  »(  l  .  (  3. 1 A 1  *i'#2T/Z2  I  •  •?♦  14 2 •  132/4. 

1.5) 

0)1?  =  0?  •  !)?•/(  TAD  1*1  2 . *C  0NS7 )  /(?.*2.»ZLTAI)l  ♦  C  0  N  S  7  ) 

UOSCM.ONS  4»»  ?»nr'L  ?•(!.♦(  ALAMDAI  I  )/(  0»C M?Cl)  ))•■?)/(.  5  « 0M2  ) 

000  II  ?  =C  (INS  4  •  (  ALAMO  A  (  1  >  /  I  0»(  .'VflO  )••?•(  X  5*  X  6  »f  ?  •  OGL  2  /flit  2 ) 

ooi)  ■  nose.*  nonir? 
go  ro  At 

37  OFl  rAC-r?».  42 /  112 

R  t  It  ’  =  2  .  •  14 2  »  r I »N S  5»KF  2 
«  e  2  =  D2  »ti 2  /  (  l  2  .  •  VS  I  AM  ) 

AL  I’ll  A2  =  AT  AO  f  (  ALAMDAI  I  )/Q  J 
CALP2  =  C0Sf  (  ALPHA2 ) 

SALP?=SINf  l  ALPHA?  ) 

C0NS9* (1 . <C At P? )*«2oFXPI  (DFLTAC ) 
CALP3'«r.ONS9-SALP2«SAlP2)/(rONS)»SALP?»SALP2) 

ALPHA  4 > AS  INFH  1  .-CALP3»CALP 1 » • • . S  > 

0U?3  =  0FLT AC • I l .♦CALP2»CAlP?/( 1 . On  1  ft *C ALP  3 »C ALP3 1 ) «C0NS3 / ( 2 . »C OSF I I 
1 ALPHA2«ALPHA3)/?. I ) 

F2*.d4/>?/Rt'H?»».2 
03* 1 .04*02 
SI42  =  14  2*02 

CONS  1 1  -  -  I  It2/r>3 

C UNS  1 2  *  (  SI4?  »C  AL  P  3/0 3  )  •  »2 
CONS  1  )»  t. 141592  7-E2/H. 

CONS  14 *S INF ( ALPHA  3) »S INF (ALPHA  3) 

0403x1  .  3 
J3*  1 

38  CONS  15* CONS  1 3«l  1 . *0403 ) / 1  040 3 ♦CONS 1 1 ) 

0403N*  (CONS  1  ?  /  I  (CONS  14/  (040  4 *0403 > -CONS  1 5 1 /  I  1.4  CONS  15  I  I  )••.?!>♦  1  . 

IF  (  AltSFI  04  03/040  3N- 1  .  )  -  .  00  1  )  4?  ,  4?  ,  3) 

39  J  3*  .1  )♦  l 

IF! 33-50)40,40,41 

40  040  4=D403N 

SOI  |F(  0403-1  .  )  4C0 ,  3  »t,  )fl 
400  040  3  =  0403  * . 1 
GO  Til  401 

41  ri  H  |  T  f  (IIIIPIIT  IAPL  42, 10  4,040  »,04I)3N,  ALPHA  » 

GO  TO  50 

42  00  3  4=1  .  5  71'7'N,  4»Ci)OS3»r  ?  •(  0?  •CALP  ,/(  11  4»CAI.P3  ))*■•’•  HI  .  *0403  )/ (l)4l)3* 
ICONS!  1  )  )  •  II  .  «  (f3?,r.ALP3/  10  4*  I  04  0  3-  I  .  )•  I  040  3-  1  .>)»••  2  ) 

Di)4'»'CnNS  J*  ( I),'  •(  Al  P?»S'32  »C  Al  P  V'  (o  I  »C  ALP  3  *03  •  (  OA  113-1  .  )  •(  0403  -  l  .  )  n  • 
1»2*(X  4  4X7»f  04  03«l. )/(Kf?**.?»(  04 1)3-1. ))> 

OODHiO?  3 ♦  i HJ  44*0045 
GO  TO  43 

4  3  f  I  HVI)?=  1  1  .  -  ODDI  -1)00)  /  (  1  .  ♦DOSE  ♦DOKI  ' 

Htllt?  ilO/l  I  M  V I )  2 
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IFI  A3SF{MTnr?/NTUT-l. )-.P2)  47,4 7, 44 

44  J2*J2*1 

I F ( J2-50)45,45,46 

45  HTOI *HTOT 2 
GO  TO  16 

46  W  X  I  T  t  OUTPUT  TAIT:  4  2  ,  11  4  ,  Ht  D  T  ,  H  T  OT  2 
GO  TO  50 

47  OOFH*X«»(  I.GF6/KE?)  ••  .2/1  1.  E3»ALAMDA(l)  »0#  13?) 

DQL*XV».G06»V2»/2»l  l  .-CONS3  I  .  5/  R2 

6T  A*  E  T  MY02/ ( 1 . ♦DOfR*OOL ) 

WHITE  OUTPUT  !  API:  4  ?  ,  1  OB  ,  RE  T  A?  ,  AMU,  AL  AMOA  (I  )  , 

1  0,02 , U2 , ZE  T A 1 , Z 1 ,  82 , RE  I , OOSF , 00 

iol.ookc.et aimp 
GO  T . I  (48,4V)  .K.ASE2 

48  WHITL  n“TPUT  f  APT  4  2 , 1 1  5  , !  C  A  SE  2  ,  RED  2 ,  RE  2 , 12  ,  DQSC ,  DQI)  l  F2 , 000 
GO  TO  491 

49  white  output  tape  42 , 1  u,  ic  ase?  ,002  3,ocj34,do45,oqd 

491  WHITE  OUTPUT  TAPE  4 2 , 1 1 7 , E T MVD? , HTO T 2 .DOCK , OOL , E TA 
IF ( TEST  )  50, 5(  ,4  V? 

492  WHITE  OUTPUT  TAPE  4?, 123 
ANNO  = ANNOO 

68  WHITE  OUTPUT  TAPE  4?, 122, ANNO 
AMO* AMOO 

69  ANPL*ANNO»AN 
AM* AMG4D0L 
\/UP*VG»AMC»ANNO 
U I P  =U I • ANNO 

CONS l  *  1 44 . »V( P/I EI*OI#0 I • AK l». 7853992) 

J*1 

Clip  3  CONS  1 

I F ( CONS l » CONS  1/44, 35- AN PSP )  /  10, 710,  700  A0545 

700  C!2P*CONSl»U.»A*ICIlP«CUP/64.35-ANPSP))*»AK  A0545 

IFI  AOSFIC  1 1P/C  I  2P-1  0001.)  710,710, 701 

701  J*J*1 

IF! J-50) 702, 70?, 703 

702  C  1 1  P * C  I  2 P 
GO  TO  700 

703  WHITE  OUTPUT  TAPE  4 2 ,  1 24 , C I  1 P , C l 2P 
C  1 1 P* ( C  1 1 P  +  C I  ?P  )  /2. 

710  OIPL*ATANr(UlP/C  1  IP  I 

DELTHI  *  AHSF(RIPl»57, 29579-81 ) 

W|P*UIP/SINHHlPl> 

CONS  l  =  C I  1 »C I  1/64.  35-ANPSP  A0545 

C0NS2=CI)P"CI IP/64. 35-ANPSP  A054S 

IF(CONSl) 704, 704, 705 

704  81=1.0 

GO  TO  706 

705  8  I  * ( 1 • ♦ A »CONS 1 1 »*AK 

706  1 F  (  CONS 2 )  70  7,  7°  7 ,  70R 

707  8  I  P  *  1 . 0 
GO  TO  709 

708  R I P* I  1 . ♦ ANCONS’ ) •• AK 

709  U2P  *U2 • ANNO 

73  CONS  l  *  I  VI I  P/II2P  )  »•? 

CONS 2’ AMO • AMO • A l AMO  A*  AL  AMOA 

CONS  3* X l PL •  ). 14! V9?  ?»C8?/(2 . «71 • I  1 .-V?) ) 

C0NS4  = 1 . - AL AMOA • AMO • TO? 
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AfllJPL  -  AMU 

J»1 

QP11-AMUPL-C0NS4 

74  ZET  A1-HC0NS2M  i.-UPll)  »»2 ) /CONS  1 1  ••  .5 
A.-1U2  -  I .  ✓  (  1 .  ♦  C  ON  S 1  •  (  l.  +  .CM/ZLTAl  )) 

QPL2-AMII2-CONS4 

I F ( AHSFlyPL 1/0PL2-I. I-.U5) 70, 7 A, 75 

75  J-JM 

I F { 3-50)76, 76,77 

76  OPl 1-0PL2 
GO  ru  74 

77  WRITE  OUTPUT  TAPE  4?  .  1 1*)  ,UPl  l ,  QPl2  ,  AMU2 
70  PSI2»Al*A2*OHTH|*A3»DF.LTiU»DELTttl 

00  INC* PS  I  2-CONS  1/1 2. -OPl 2) 

WAVtU22.5-IC(JNSl*C0N$2M  l.-0Pl?)-»2) 

I F (  ICASF-4I0C, 01,00 
GO  PS  IC*o.O 
OOCAV-O.O 
GO  TO  R4 

81  DoCAV  2  PSICMR IP-l. )-C0NSl/UPL2 
84  WAVE=WAVEU2--.S-U2P 

«Ei=OHYm»wAve/i  i2.»vstapl) 

FI- .0462/KE  I  «  •  .  2 

0USF  =  X2PL  -F 1 -0FLl-WAVrU2  / (OMYDi -QPL2 ) 

WIU2*CONSl--.5 

CL-COSFIUIPl> 

DeiT«l.-ZETAl*0Pl?-ICl*CH2)/l2.-WlU2*IZl-ll.-V2>/5.14159274?.-V2)> 
DOOl2X3PL-OFLT-DElT  -C0NS1/QPL2 
A.'*  AMO 

00RC*X4PL •OCL  T  / ( ALAMOA- AM) 

000  IN- PS  I 3-1 ( UPL2-0I --2-I Al  AMOA- AM- ALAMOA)  ••  ?  )  /  1 2 . -QPI.2  ) 

DQOPL  =000- AM -Q-  ( ( C0NS2  +  0PL2 -QPl  2 )  *0/ II  ALAMOA-  AL  AMDA  ♦Q»0  )  -QPL2  ) )  --. 
15/QPL2 

HTOT-QPL2-U2P-U2P/32. 175 

ETAH-I  1 . -DO  I  NC- OOC  A  V-OQUL  -  001)  I  N-OQDPL  )  /  1 1 . -DQSF  *0080 

HO-ETAH-HTOT 

DELTP-GAMP-HO/I 44. 

RE2-U2P-D2/I  12.* VST  API) 

DyFR*X8PL-( 1 . E6 /RF  2 ) •* • 2/ I  1 . E 3- ALAMOA* AM0-GPL2-B2 ) 

DQL Pl-X^Plj  -  ,0Q6-V?«v2-(  I2.*QPL2-CONS2*OPL2*QPL2) / 12. -QPL2) )  --.5/B2 
ETA-ETAM/  I  1 .  ♦OOFJUOOIPI  ) 

HP-HO-AMO-VO-GAMP-ANNO/ IFTA-550.  ) 

ANS-AMO-ANNO 

WRITE  OUTPUT  TAPE  4 0 , 12 0 , VO , OEL T P , E T A , ANS , ANNO 

WRITE  OUTPUT  TAPE  42, 121, AMO  , OPL 2 , OF l TB I , F T AH, MO ,0C LTP , E T A , HP 
IFIAMQ-AMN  (06,07,87 
06  AMO*AMO»OELTM 
GO  TO  69 

87  IF  I ANNO-A NNON ) fl 8 , 8 9 , 89 

88  ANNO*  ANNO  *OEl T  AN 
GO  TO  6fl 

89  CONTINUF 
50  CONTINUF 

GO  TO  l 
ENO 
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